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Abstract 

Asymptotic-in-time feedback control of a panel interacting with an inviscid, subsonic flow is considered. 

The classical model [21] is given by a clamped nonlinear plate strongly coupled to a convected wave equation 
on the half space. In the absence of imposed energy dissipation the plate dynamics converge to a compact 
and hnite dimensional set | 16l I17 |. With a sufficiently large velocity feedback control on the structure we 
show that the full flow-plate system exhibits strong convergence to the stationary set in the natural energy 
topology. To accomplish this task, a novel decomposition of the nonlinear plate dynamics is utilized: 
a smooth component (globally bounded in a higher topology) and a uniformly exponentially decaying 
component. 

Our result implies that flutter (a periodic or chaotic end behavior) can be eliminated (in subsonic flows) 
with sufficient frictional damping in the structure. While such a result has been proved in the past for 
regularized plate models (with rotational inertia terms or thermal considerations |1411321 [3^ I40p . this is 
the first treatment which does not incorporate smoothing effects for the structure. 

Key Terms: strong stability, nonlinear plates, flow-structure interaction, PDE with delay, attractors 
2010 AMS: 74F10, 74K20, 76G25, 35B40, 35G25, 37L15 

1 Introduction 

This work is motivated by physical problems in aeroelasticity. Flutter of a thin structure immersed in a flow 
field occurs (as a bifurcation at a particular flow velocity) when the structural modes couple with the dynamic 
loading of the flow, creating a feedback 0 m]. The canonical example of panel flutter occurs in aircraft 
paneling: panel deformations perturb the flow field, which then determines the dynamic loading on the plate 
(resulting in structural deformation, etc.). The principal challenge in aeroelasticity is the control, suppression, 
or prevention of the flutter of mechanical structures in flight [23] . Though the flutter point may be 

predicted from linear dynamics, accurately describing post-flutter displacements of a panel requires a nonlinear 
model pn 122] . 

From a mathematical point of view the flutter problem may be described as stabilization to a given station¬ 
ary set for a coupled evolutionary PDE system: a nonlinear clamped plate interacting with a potential flow 
0I1I1I1H]- Here, we consider Berger plates in the absence of rotational inertia (consistent with the conventions 
of aeroelasticity 0I1IJI11I13] and references therein), where the plate is thin and in-plane accelerations are not 
accounted for [3D]. The task at hand, then, is to determine conditions under which the PDE system is strongly 
stable, in the sense that full flow-plate trajectories converge (in the strong topology of the underlying finite 
energy space) to a set of interest (e.g., the equilibria set for the dynamical system generated by the solutions). 
Such a strong convergence to equilibria does not seem possible in the absence of dissipative mechanisms. Thus, 
we seek an interior velocity feedback control (of sufficient size—depending on the intrinsic parameters of the 
problem, not on the particular initial data). Under some physical assumptions this will imply that flutter is 
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Air Force Office of Scientific Research with grant AFOSR-FA99550-9-1-0459. J.T. Webster was partially supported by National 
Science Foundation with grant NSF-DMS-1504697. 

f University of Memphis, Memphis, lasiecka@memphis.edu 
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non-present asymptotically in time. Our result demonstrates that flutter can he eliminated (in the subsonic 
regime), provided a sufficiently large damping is present in the plate. In agreement with the physics of the 
problem, the linear model for the flow dynamics remains undamped. 

In this treatment, we primarily focus on the Berger plate model. Other work on flutter dynamics have 
considered Berger-like beam and plate models (see |13j and references therein), as well as the scalar von 
Karman equations [H [la na mi- Many of our supporting results here are valid for both the von Karman 
and Berger plate, but the principal results presented below in Section [1.61 concern Berger plate dynamics. The 
Berger plate is a nonlocal, cubic-type nonlinearity. The Berger nonlinearity is arrived at via a simplification of 
the scalar von Karman equation 1351143] (with the scalar von Karman equation itself a simplification of the 
full von Karman system |291 130j lrl However, there are times (such as in the main result presented here) when 
working in higher topologies for the dynamics produces a stark difference in the behavior of the von Karman 
plate and the Berger plate. 

Physically, when the plate has no portion of its boundary free, the Berger approximation is taken to be 
valid [23 132 ] (and especially in the case of clamped conditions, as we take here). The von Karman plate 
dynamics were considered in (22] (and in the treatments leading up to [22]! including [121 |39l [40] and [14] 1. 
In |32] . large static and viscous damping was considered to obtain the main result (which is valid for both 
the non-rotational von Karman and Berger dynamics). Here, we consider the more physical case where only 
viscous damping is active in the plate dynamics (though we cannot dispense with the size requirement). 

1.1 Previous Literature in Relation to the Present Work 

The study in m provides long-time asymptotic properties of (finite energy) solutions for the model herein 
(for Berger, von Karman, and other nonlinearities of physical type |18] 1 with no restrictions on the damping 
size —it can be taken to be zero. In this case, the originally “rough” plate dynamics become asymptotically 
regular and finite dimensional in nature. Thus, the flow has the ability to dissipate mechanical plate energy 
and thereby induce a degree of stability to the structural dynamics. However, this stability, in general, may be 
of a chaotic character. It is one of the principal purposes of this treatment to show that, in the subsonic case, 
sufficiently large dissipation in the structure alone prevents this chaotic behavior and leads to strong stability 
to the equilibria set. 

We now mention three other closely related scenarios in the literature to date: (i) the addition of ther¬ 
moelastic dynamics to the nonlinear plate, (ii) the presence of rotational inertia terms and strong mechanical 
damping, and (iii) large static and viscous damping in the model presented here. (These treatments address 
von Karman plates, though the results also hold for Berger plates.) The treatments in [ 211110 ] consider the 
plate with the addition of a heat equation in the plate dynamics. In this case no damping is needed, as 
the analytic smoothing and stability properties of thermoelasticity provide ultimate compactness of the plate 
dynamics and, furthermore, provide a stabilizing effect for the flow to equilibria (in subsonic flows). As for 
(ii) mentioned above, results on ultimate compactness of plate dynamics, as well as convergence (for subsonic 
flows) of full flow-plate trajectories to points of equilibrium, were announced in |12] . with a proof appearing 
in [14]. The recent treatment [21] addresses (iii); it is shown that by considering “large” static damping and 
“large” viscous damping—of the form D ■ [ut + u] in the plate equation, D sufficiently large—flow-plate tra¬ 
jectories with finite energy initial data converge to the set of equilibria. The techniques in [32] are rooted in 
Lyapunov methods and utilize critically the prior work in [miniEH ED]. The key point here is that for large 
static and viscous damping, exponential decay of plate velocities provides global-in-tiine Hadmard continuity 
which permits approximation by smooth data (for which the desired result holds—Theorem 14.91) . 

The main goal of the present work is to obtain strong stability results ( strong convergence to a set of 
equilibria) without accounting for smoothing effects in plate dynamics and without without assuming large static 
damping imposed on the structure (plate). The primary mathematical issue to contend with is low regularity 
of the hyperbolic Neumann map (boundary into the interior)—i.e., the failure of the uniform Lopatinski 
conditions in dimensions higher than one. This precludes direct analysis of the coupling via the compactness 
properties of the aforementioned Neumann map. (Such an approach was used critically in previous analyses 
[H [El EH [H EH ED].) In fact, while we still rely on these past developments, the key point of departure is that 
we can no longer afford to treat the problem component-wise. Rather, we must rely on global analysis involving 
relaxed long time compactness [201141] , which depends on uniform long time invariance of higher energies [32] . 

^In fact, in many cases the analysis of the scalar von Karman plate subsumes that of Berger m 
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This latter property is highly non-trivial due to the effects of the plate nonlinearity combined with the lack 
of damping exhibited by the flow. In order to overcome this obstacle we shall rely on the idea introduced in 
[13] and used in m where structural equation under the influence of a flow is represented (for large times) 
by delay system. It is then shown in m that the corresponding decomposed part admits a smooth global 
attracting set. In order to extract a long time convergence valid for the entire flow-structure system it would 
suffice to have a control of Li rate of convergence of delayed plate solutions to such attractor. This brings 
forward a concept of exponential attractor, which by definition provides exponential convergence rates, but it 
may loose ’’smoothness ” properties. The above leads to a well known dichotomy between smoothness and 
rates of convergence to attracting sets. Resolving this dichotomy lies in the heart of the problem presented. A 
delicate argument based on a careful analysis of delayed dynamics along with an uniform Hadamard continuity 
property exhibited by of a suitable decomposed part of the full dynamics only, allows to reconcile this dilemma. 
As a consequence, we show that, with a sufficiently large viscous damping parameter, the entire flow-plate 
finite energy dynamics exhibits strong convergence properties. We emphasize that our results require only 
viscous damping in the structure and do not make use of any advantageous parabolic (smoothing) effects or 
additional static damping. One may argue with the requirement of sufficiently large size of the damping (not 
needed in parabolic like problems), however the fact that we de facto stabilize delayed structure (with inherent 
instability) justifies mathematically |36l ES] a need for large viscous damping. 

Our result provides mathematical insight to the panel flutter problem; in |16j the flow, while driving the 
plate dynamics, also contributes a stabilizing effect to the plate dynamics. Our present result indicates that 
for a panel in a subsonic flow, strong stability of the plate (due to mechanical damping) can be transferred 
to the flow, in some sense. (We note that in all real panels there is some degree of internal damping.) These 
results are in agreement with experimental and numerical studies wherein divergence (or ‘buckling’) of panels 
is observed for subsonic flows, in contrast to chaotic or periodic end behavior (z.e., flutter) in the case of 
supersonic flow velocities. We quote from a recent survey article authored by foremost aeroelastician E. 
Dowell [24] with regard to subsonic flows: “...if the trailing edge, as well as the leading edge, of the panel is 
fixed then divergence (static aeroelastic instability) will occur rather than flutter.” In the analysis below, the 
subsonic nature of the flow is critical, as it provides a viable energy balance equation; this energy balance is 
not available for supersonic flows [18]. 

The engineering literature in dEillS] analyzes the validity of Berger’s plate approximation (originally 
appearing in 0) and concludes that it is accurate in the case of clamped and hinged boundary conditions. 
The approximation is based on the assumption that the second stain invariant is negligible. Mathematically, 
the Berger dynamics are addressed in (TOj US] |T9] (in the context of aeroelasticity), and also in the papers 
[261 EU- Much of the abstract work in [14] (which is focused on the scalar von Karman nonlinearity) also 
applies to the Berger dynamics. The main results appearing in [m nil [mill] for the model presented above 
discuss and analyze both the Berger and von Karman nonlinearities, as well as other “physical” nonlinearities 
satisfying a specific set of bounds which are referred to as nonlinearities of physical type). 

1.2 Outline 

The remainder of the paper is organized as follows: 

[§1] The remainder of the Introduction is devoted to a discussion of the mathematical model, including the 
relationship between our analysis and previous literature. We discuss the fundamental notions of well-posedness 
and associated energies, and describe some key results from dSlEl] which motivate this analysis. We then 
describe the main results. 

[§2] Section [2] describes the technical preliminaries needed to read the proof of the main results. This section 
can also be read as a detailed summary of the preceding work on the specific ffow-plate model described herein. 
[§3] Section [3] describes the supporting results from [16], used here, which deal with the existence of a compact 
global attractor for the plate dynamics, as well as some weak stability properties. 

[§4] In Section [4.2l we prove the primary supporting result upon which our later analysis is based: Theorem 14.31 
which says that a global-in-time bound on the plate dynamics in higher norms yields the desired convergence 
to equilibria result for the full flow-plate dynamics. This section outlines a strategy of passing convergence 
results from smooth initial data onto finite energy initial data, and we describe in Section 14.51 the primary 
technical obstacles in showing the convergence to equilibria result for finite energy initial data in the absence 
of smooth effects or large static damping. 
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[§5] Section [5] begins by providing an overview of the main steps in the proof, and their relation to the 
established results in Section [21 (We believe this to be helpful to the reader, owing to the technical, and 
multistep nature of the proof.) The proof of the main result follow in substeps: a decomposition of plate 
solutions (Sections 15.31 and 15.4L a key decomposition for the full flow-plate dynamics, an analysis of the 
decomposed dynamics, and a concluding argument which yields the final result. 

[§6] In Section [5] we discuss the open problem of eliminating the size requirements on the damping. 

[§7] Readers will find the authors’ acknowledgements in the final section. 


1.3 Notation 

For the remainder of the text we write x for {x,y,z) G or (x,y) G ft C as dictated by context. 

Norms || • || are taken to be L 2 [D) for a domain D. The symbols u and r will be used to denote the unit normal 
and tangent vectors to a given domain, again, dictated by context. Inner products in L 2 (R+) are written (•, •), 
while inner products in = 9R+) are written (•,•). Also, H^{D) will denote the Sobolev space of order s, 

defined on a domain D, with Hq{D) denoting the closure of C^{D) in the H^{D) norm (which we denote by 
II ■ llff®(D) or II • ||s__d). When s = 0 we may abbreviate the notation to || • ||. We will also utilize the homogenous 
Sobolev space 

iiwii.V^j. 

We make use of the standard notation for the trace of functions defined on i.e., for (f) G il^(]R|^), tr[(j)] = 
restriction of cj) on the plane {x : z = 0}. (We use analogous notation for tr[w\ as the map from 
H^{n) to Finally, we will denote a ball of radius R centered at c in a topological vector space X 

by Bx{c,R); when the center is the zero element we simply write Bx{R)- 


{ k-i 
(/.(x) € 4-(R^) : ml,, = ^ 

j=o 


1.4 Mathematical Model and Natural Energies 

The gas-flow environment is modeled by = {{x,y,z) : z > 0}. The plate is immersed in an inviscid 

potential flow (over body) with velocity C/ < 1 in the x-direction. (Here we (i) normalize 1/ = 1 to be Mach 1, 
i.e., 0 < C/ < 1 corresponds to a subsonic flow, and (ii) consider only the dynamics in the upper half space.) 
The panel is modeled by a bounded domain H C = {x : 2 ; = 0} with smooth boundary dft = F, and the 

scalar function u : fix R+ —>■ K. represents the (transverse) displacement of the plate in the z-direction at {x, y) 
at the moment t. For the flow component of the model we make use of linear potential theory [6l|23], and 
the (perturbed) flow potential is given by ()) : M3 X IR+ ^ R. The strong coupling occurs (i) in the dynamic 
pressure term (RHS of the plate equation), which contains the acceleration potential of the flow, and (ii) in 
the downwash of the flow (Neumann condition), the latter including the material derivative of the structure. 


'utt + -I- kut + f{u) = po + rQ,tr[[dt + Uda)4>] 

U(0) = Mo; Mt(0)=Mi, 
u = dvU = 0 
{dt + Udj:f4> = A4> 

<^(o) = 0o; </>t(0) = (/>i, 

[{dt + Udm]^^t 

The notation xq corresponds to the restriction (of the function’s domain) to H (in the appropriate functional 
sense) for functions supported on aR3 ; correspondingly, we have the extension by zero for functions supported 
on n to the entirely of (in particular, for functions in Uq (H)), denoted by Uext- The term po is a static 
pressure on on the plate. The nonlinearity of principal interest here is the Berger nonlinearity: 


in n X (0,T), 

on F X (0,T), 
in X (0, T), 

on X (0^^)- 


/b(m) = [6-||Vm|P]Am. (1.2) 

The parameter 6 > 0 is a physical parameter [35l [43] which corresponds to in-plane bending or stretching!! 

■^In full generality, b can be any real number; we choose the more mathematically interesting—non-dissipative—case, with 
fc > 0. 


4 




(1.3) 


The plate energy is defined as usual [I1I30]: 

E,,iu)=l[\\u,\\l + \\Au\\l]+U{u). 

n(u) is a potential of the nonlinear and nonconservative forces, given by 

Il{u) = nB{u) = i||Vu||^ - ^||Vu||^ - {po,u)q. (1.4) 

The natural energies associated with the subsonic wave and interactive dynamics are given below: 

Efiicj,) = (1.5) 

Eintiu,(l)) = 2U{tr[4>\,Ux)n- (1.6) 

The total (unsigned) energy is then defined to be 

£{u(t),(l){t)) = £{t) = Epi{u{t)) + Efi{4>(t)) + E,„tiu(t),4){t)). (1.7) 

We will also need to consider positive energies, so we define 

^4u) = i||Vu||^ i5.(u) = i[|K||2 + ||Au||2] + n4u), (i.s) 

£^{u,4>) =E^{u) + Efi{4>). (1.9) 

According to these norms, the natural energy space is then: 

Y = YpiX Yfi = {Hl{n) X L2{n)) X (w 1(]R3 ) x L2(R+)) , (I-IO) 

taken with norm 

||(u,z;;(/),V’)||v = ||Au||^ + ||u||^ + ||V ((>||23 + HV'IIr^- (1-11) 

We will also consider a stronger space: 

W = {Hl{n) X L2{n)) X X L2(M^)) . (1.12) 


1.4.1 Damping 


The damping coefficient is given by fc > 0, representing viscous damping (or dissipation) in the model. The 
can be viewed as intrinsic damping for the structure, or imposed as a velocity feedback control. For certain 
results (see |16j l we can consider fc = 0, however, for the main results below we will need specific assumptions 
on the size of the damping coefficient k. 

Additionally, in some auxiliary results presented below (especially those related to [32]) we will refer to 
static damping —this occurs when the viscous damping kut is accompanied by a term of the form j3u. The 
resulting system, which will be referred to below, is 


Utt + A'^u+ kut + Pu + f{u) =po+rQtr[{dt + Udx 


u{0) = Uq] Ut{0) = Ui, 


u = di,u = 0 

{dt + UdxYP = A(j) 

(()(0) = 0o; (/>t(0) = </>!, 

d^(l)= [{dt + Udx)u]^^^ 


in n X (o,r), 

on r X (0, T), 
in X (0, T), 

on ®(x.y) X (0>^)- 


(1.13) 


Note that, when /? = 0 this system collapses to (HD), the principal model under consideration here. In this 
case, we adapt the plate energy Epi and augment it by the quantity /3||it|p: 


Epi,f!{u) = i [llutll^ + ^||m|P + ||Am|P] +n(u) = Epi{u) + ^WuW^. 

In fact, the primary results in [32] consider large static and viscous damping for the von Karman plate. The 
main aim of the present paper is to obtain strong stability result without any static damping, i.e. with /3 = 0. 
This will be achieved by asserting that sufficient amount of stability can be harvested from the flow alone. 
However, rigorous proof of this physically plausible fact is quiet involved and delicate. 
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1.5 Well-posedness and Fundamental Notions 

The dynamics discussed above, when cast in the appropriate framework, are well-posed 

(For precise definitions of strong, mild, and weak solutions for the flow-plate dynamics consult [T^[T^l32ll^ .l 
The following specific well-posedness result is established in [HIH]: 

Theorem 1.1 (Nonlinear Semigroup). Assume U < 1, po G ^ 2 (^ 1 ). Then for any T > 0 (11.11) has 
a unique strong (resp. generalized — and hence weak) solution on [0,T], denoted by Stijjo), for initial data 
j/o = {uq,ui;4io,4>i) G Y. (In the case of strong solutions, the natural compatibility condition must be in force 
on the data 

dz4>Q = (ui + Uuox)ext-) Moreover, {St,Y) is a (nonlinear) dynamical system. Additionally, for j/o S Yg, 
{St,Ys) is also a dynamical system. Weak (and hence generalized and strong) solutions satisfy the following 
energy equality. 

£{t) + k f [ \ut{T)\^dTtdT = £{s) (1-14) 

J s Jq 

fort > s. Moreover, the solution St{yo) is stable in the norm ofY, i.e., there exists a constant C'(llyollv) such 
that for all t > 0 we have: 

||^t(yo)l|y <C'(||yo||F). (1.15) 

In addition, the semigroup St is locally Lipschitz on Y 

\\St{yi)-St{y 2 )\\Y <C{R,T)\\y^-y 2 \\Y. ^MIy < R, t<T (1.16) 

Remark 1.1. The proof of Theorem ll.ll is based on a suitable rescaled variant of monotone operator theory and 
is given in [44]. See also m for a different approach based on viscosity method. In the case when U > 1 the 
situation is more delicate and dealt with in [TB] . In this latter case (11.141) no longer holds and the inequality in 
(11.151) is valid only locally, for t < T . This is the result of the loss of dissipativity occurring in the supersonic 
case. 

For the above semigroup we introduce the dynamics operator T : ^(T) C Yg For its precise structure, 

we give reference to [T5| |TH| [44] which provide the details of the abstract modelH A key property needed in 
this treatment is: 

^(T) c (iJ'^ (h i?o)(ll) X H^{n) X H^(Rl) X H\Rl). (1.17) 

Remark 1.2. As per [Ml [151144] th® natural invariance of the dynamics is with respect to the norm || •||y. 
However, via the hyperbolic-type estimate 

ll<(’(l)llL2(Ri).) - ll'(^o||L2(Ri).) + y \\4>t{T)\\L^(^R3^)dT, (1.18) 

invariance in Yg can be recovered on finite time intervals. 

In order to describe the dynamics of the flow in the context of long-time behavior it is necessary to introduce 
local spaces, denoted by Yfi^pi 

11(00, <('i)||yj.,.p = f |V0o|^ + |0i|^dx, 

JKg 

where ATp = {x € Ri].; |x| < p}. We denote by Yp the space Ypi x Yfi^p, and we will consider convergence (in 
time) in Yp for any p > 0. Such convergence will be denoted by convergence in the sense ofY. 

By virtue of the Hardy inequality [HI p.301] 

\\4>o\\l2{Kp) ^ C'p||V0o|li2(R3.) 


and hence 

11(00, 0l)||v^i,p < 11(00,0l)llff ^(Kg)xL2{Kg) ^ 11(00, 0i)||y^,. 

^We suffice to say that Ball’s method provides the generator of the nonlinear semigroup with appropriate dense domain ^(T). 


6 







Remark 1.3. We clarify the relation between the various flow topologies: H Yfi, and Y/i,p. 

Initial flow data {(po, <j)i) are typically chosen in x L 2 (K+)- Solutions (u(t), Ut(t); <^(1), are global- 

in-time bounded in the topology of Y, but not necessarily global-in-time bounded in the full Yg norm (owing 
to the contribution of the L 2 norm of <('(!))• Hence, for considerations involving t —^ 00 , we require a sptially- 
localized flow perspective. By restricting to any ball Kp C for the flow dynamics we obtain global-in-time 
boundedness in the topology of Y as well as local compactness results. 

By the boundedness in Theorem ll.il the topology corresponding to Yp (for any p > 0) becomes a viable 
measure of long-time behavior which we will refer to as the local energy sense, and it is appropriate to take 
limits of the form lim ||S't(j/o)||F„, for any p > 0, where || • Hv = || • and St is the flow on Y 

i—>-oo ^ ^ P J 

associated to the well-posedness in Theorem II.II above. 

We now note the boundedness (from below) of the nonlinear energy. Such a bound is necessary to obtain 
the boundedness of the semigroup quoted in Theorem 11.11 above. First, we have [HI Lemma 5.2, p. 3136]: 

Lemma 1.2. Let the hypotheses of Theorem 11.11 be in force. Then for generalized solutions to HH), there 
exist positive constants c, C, and M positive such that 


cS,,{t) - Mpg^b < £[t) < C'f,(t) -b (1-19) 

The proof of Lemma 11.21 given in [44] relies on two estimates contraolling lower frequencies. Since these 
estimates will be used frequently in the sequel, we shall formulate the relevant results. The first estimate 
controls interactive energy Eint on the strength of Hardy inequality (see [B p. 301]), 

Lemma 1.3. 

\E^nt[t)\<U5\\ym\\li + ^\\u.ml. ^> 0 , (1.20) 

The second one controls lower frequencies [Ill p. 49] by exploiting nonlinear effects of potential energy: 
Lemma 1.4. For any u € H^(fl) fl ilg (H) and e > 0 there exists such that 

||u 1 [q < e[llA u 1 ]q + n*(u)] + Mg. 

Remark 1.4. In fact, in the case of Berger’s nonlinearity Lemma [L4l holds without the term jjAttjjQ. However, 
in the case of other nonlinearities, such as von Karman, the presence of A is indispensable. For this reason we 
cite a more general version of this lemma which can be used for von Karman models as well. 

Note: global-in-time boundedness of solutions cannot be obtained without accounting for nonlinear effects. 
From the above lemmata and energy inequality we also have (see [151 l44] l: 

Lemma 1.5. Let the hypotheses of Theorem \l.l\ be in force. Then any generalized (and hence weak) solution 
to (dl) satisfies the bound 

sup|llutl|^-b llAu]]^ -b l|(/>t||R3 + ||V()-1|L I < C'(1]2 /o||y) < + 00 . (1.21) 

#>n + + J 


Moreover, from energy inequality 11.141 and (11.211) the following dissipation integral is finite (see [B1321 1: 

Corollary 1.6. Let the hypotheses of Theorem ] 1.1\ be in force and assume k > 0; then the dissipation integral 
is finite. Namely, for a generalized solution to (dl), we have 



<K < 00 , 


where K depends on the particular trajectory. 

We now state and discuss the stationary problem associated to Id, which has the form: 


A^u-b /(u) = po(x) -b UrQtr[da;4>\ 

X € H 

e 

II 

e 

II 

0 

xe F 

A(j) - U^dlf = 0 

X e 

— U QxUe-Kt 

X G 


( 1 . 22 ) 


7 







This problem has been studied in long-time behavior considerations for flow-plate interactions, most recently 
in [141 Section 6.5.6]; in this reference, the following theorem is shown for subsonic flows (this is given as [141 
Theorem 6.5.10]): 


Theorem 1.7. Suppose 0 < U < 1, k > 0, with po G L 2 {^). Then weak solutions (tt(x), (/)(x)) to (11.2211 exist 
and satisfy the additional regularity property 

{u, (f) G (d Hl){n) X ). 

Moreover, the stationary solutions mentioned above correspond to the extreme points of the potential energy 
functional 

1 1 t/^ 

D{u,(i)) = -||Au||^-bn(M) -h -||V((i|]2^ - —-bC/(5,;U,fr[(/)])a, 
considered for (u, </>) G Hq{Q) x IT^(K^). 

The potential energy D{u, (f) is smooth on the space x and moreover 


D{u,cf)>c ||Au||2 +||V</>|lL 


-C. 


This latter property is a consequence of Lemma |1.41 Thus, it achieves its minimum, and the extremal set of 
the functional D is non-empty. We denote the set of all stationary solutions (weak solutions to (11.221) above) 
as M, that is 

M = {{u,^) G Hq{TI) X lT^(]Ri]_) : (u,(^) satisfy (11.221) variationally}. 


1.6 Statement of Main Result 

Theorem 1.8. [Principal Result] 

Let 0 < [/ < 1 and assume po G ^ 2 (^ 1 ) and b > 0. Assume yo = {uq,Ui', (po, (pi) G Y. Then there is a 
minimal damping coefficient kmin > 0 (not depending on the particular solution) so that for k > kmin > 0 any 
generalized solution (u(t),p(t)) to the system with localized initial flow data (i.e., the supports of (po and pi 
are contained in some Kp„) has the property that 

= 0, for any p > 0. 

Remark 1.5. In this result, and all results below, the minimal damping coefficient kmin depends on the loading 
Po and b, as well as the domain LI, the flow support parameter po, and the unperturbed flow velocity U, but 
are independent on the particular initial data of the system. 

If we make a further physical assumption that Af is an isolated set (e.g., finite), we have the following 
second main result as a corollary: 

Corollary 1.9. Assume that Af is an isolated set. Let the hypotheses of Theorem \l.S\ be in force; then for any 
generalized solution (u,p) to (11.11) (with localized flow data, as above), there exists a solution (u,p) satisfying 
(|1.22l) such that 


{lk(0 “ + ||Mt(i)|li2(n) + ||<(>(i) - + ll?^t(i)||i2(Kp)| 

= 0, for any p > 0. 

Remark 1.6. For given loads b and po, the set of solutions is generically finite. This is to say that there is an 
open dense set TZ C L 2 {Lt) x R such that if {po,b) G TZ then the corresponding set of stationary solutions Af 
is finite. This follows from the Sard-Smale theorem, as shown in m Theorem 1.5.7 and Remark 6.5.11]. 

In addition to the primary results, we mention some results obtained along the way which are of interest in 
their own right. Owing to the necessary background information to precisely state these results, at this stage 
we describe them only informally. 












First, in Section 14.21 (following a few technical descriptions and lemmas) we prove the main supporting 
result: Theorem l4.3l This result says that global-in-time bounds on plate dynamics in higher norms yields the 
desired convergence of the flow-plate dynamics at the finite-energy level (in Yp, for all p > 0). It appeared in 
|32j as a step in the proof of the main result, but was not presented independently. 

In Section 1131 we state and prove Theorem 14.71 This result shows that for “smooth” initial data for the full 
flow-plate dynamics we can propagate regularity on the finite time horizon. This result (or slightly modihed 
versions of it) will be used many times below. For the plate dynamics alone 

this regularity can be propagated on the infinite time horizon (with corresponding uniform bound). More¬ 
over, this result does not depend on the size of the damping. 


2 Technical Preliminaries 

In this section we discuss a few key technical results and constructions which will be used critically later. 
References are provided for each of the considerations below. 


2.1 Flow Potentials with Given Neumann Data 


In what follows it will be necessary to consider the flow equation with given Neumann data, 
problem: 


{dt -I- Udxf(j) = A4> 


d^4> 


z=0 


= /i(x, t) 


in 

in 




— 4>0] — 4’! 


Consider the 


( 2 . 1 ) 


We have the following theorem from [3 [13 [34]: 


Theorem 2.1. Assume U > 0, U 1; take {(l)o,4i) S x L 2 (IR.^). If 

h£ C ([to, oo); il^^^(K.^)) then (12.111 is well-posed (in the weak sense) with 


4 G C([to,oo);iJi(R3 )) , G C([to,oo);L2(Kt)) • 

Remark 2.1. In fact, a stronger regularity result is available. Finite energy x ^2(^1)) solutions are 

obtained with h G il^/^((0,T) x R^) [STJESj- However, the corresponding estimate does not have sufficiently 
good control of the dependence on time T, and hence is of limited applicability in the present (stabilization) 
context. 


2.2 Flow Superposition and Compactness Criteria 

We may decompose the flow problem from into two pieces corresponding to zero Neumann data, and 
zero initial data, respectively: 


where we will later take: 


(dt + Ud,)^4* = Ar 

4*{0) = 4o; ,/)t(0)=())i 

A4* = o 


in X (0, T) 
in aiRi). X (0, T) 


{dt + Ud^)'^4** = in Rl x (0, T) 

4 **{ 0 )= 0 ] ( j )(*( 0)=0 

du4** = h{x,t) in x (0,T) 


h{x,t) = -[Wt -b C/U3,]ext G C'(T2(R^)) 


( 2 . 2 ) 


(2.3) 


(2.4) 


Following the analyses in [31 [13 [33 [40] we denote the solution to (12.21) as 4* and the solution to (ESI) as 
cj)**] then, the full flow solution cj) coming from (EU has the form 


4{t) = 4*{t) + 4**{t)- 
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Remark 2.2. The analysis of (j)* below is identical to that given in [3 mm [35] ■ However, the treatment of the 
second part which corresponds to the hyperbolic Neumann map, is very different than these references 

due to the loss of a derivative in the Neumann map. Indeed, with rotational inertia in place (ut € H^(fl)) one 
has (for finite energy solutions) h € C(0,T; On the other hand from [34] h G i—>• (/>** G 

n C'^(L2(R+))—where the latter is of finite energy. Thus the Neumann map is compact in this 
case. In the absence of rotational inertia one has only that h G (7(^2 (R^))- Such boundary regularity does not 
reproduce finite energy solutions, with a maximal regularity (j)** € fl yielding 

the loss of 1/3 derivative. This loss is optimal and cannot be improved [13] • This fact underscores that the 
entirely component-wise analysis of finite energy solutions for the flow-plate model successfully performed in 
the past literature [111 EH EO] cannot be utilized here. 

For the analysis of (j)* we use the tools developed in mm- Using the Kirchhoff type representation for the 
solution t) in (see, e.g.. Theorem 6.6.12 in [l4])i can conclude that if the initial data i/o and (pi are 
localized in the ball iFp, then by finite dependence on the domain of the signal in three dimensions (Huygen’s 
principle), one obtains </*(x, t) = 0 for all x G Kp and t > tp. Thus we have that 

{dt -G Udj;)tr[p*] =0, xgO, t>tp. 

Thus p* tends to zero in the sense of the local flow energy, i.e., 

(t)\\l2{Kp) + \\4'ti't)\\L2{Kp) ^ t — >■ OO, (2.5) 

for all fixed p > 0. 

We now introduce a compactness criterion in the local flow energy sense. 

Lemma 2.2. Let {{p^,PT)}m ® bounded sequence in Yfi = IUi(]Ri^) x L 2 (K+) and let rj > 0. If for any 

p > 0 there exists an N{p) and C{p) so that 

\W<^e\\lKp + UT\\lKp<C{p) 

for all m > N{p) then the sequence , PT)}m compact in Yfi. 

This is given as Lemma 10 in [39l p. 472] (where it is proved) and is utilized in a critical way in [40| as well. 

For the term p** we have the following theorem which provides us with an explicit form of the solution 
(for a proof, see for instance [T4[ Theorem 6.6.10]). Below, we utilize the notation: 

p{x,t,s,9) = f{x- Ki{0,s,z),y- K2{9,s,z),t- s ), 


and 

Ki(9, s, z) = Us + 'J s’^ — z"^ sin0, K 2 ( 6 >, s, z) = \J s"^ — z"^ cos 9. 

Theorem 2.3. Considering the problem in (|2.3I) with zero initial flow data, and considering 

/i(x,t) = -[ut{x,y,t) + Uux{x,y,t)\exu 

there exists a time t*{Ll,U) such that we have the following representation for the weak solution p**{f) for 
t>t*: ^ 

p**{x,t) = - / / {ul{x,t,s,9)+Uul{x,t,s,9))d9ds. (2.6) 

where x(s) is the Heaviside function. The time t* is given by: 

t* = inf{t : x(t/, 0, s) p. ft for all {x,y) G II, 9 G [0, 27 r], and s > t} 

with x([/, 9, s) = {x — (U + sin0)s, y — s cos 9) C (not to be confused with x having no dependencies noted, 
which is simply x = {x, y)). 

Moreover, we have the following point-wise formula for the derivative in t m p. 4^0](which is justified 
for plate solutions with initial data in ^(T), and can be taken distributionally for data in Y): 


1 1*277 1*277 

'/r(x,<) = ^{y ul(pi,t,t*,9)d9 - J u\{yL,t,z,9)d9 


(2.7) 


Jo 

with Mg = sin 9dx + cos Ody. 


[dxul] (x, t, s, 6)d6ds + 


fz JO 




:[MeuJ](x, t, s, 9)d9ds'^ 
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We note that with (0o,<^i) S x L2(M+) one obtains [T5115i] 

{(j)*(t), (t)) € X L2(M+). Thus, by Theorem ll.il we also have that 

(<^**(t),</.r(t))eTf'(Kt)xL2(K^). 

Remark 2.3. Note that this last property is not valid for a flow solution with h € i2(R^) boundary Neumann 
data. The general theory will provide at most 

X [0,T]) [42]. However, the improved regularity is due to the interaction with the plate and the 
resulting cancellations on the interface. Moreover, we also obtain a meaningful “hidden trace regularity” for 
the aeroelastic potential on the boundary of the structure [niiiij: 

{dt + Ud,)tr[cj)] € L2(0, T; (2.8) 


where T is arbitrary. 

Additionally, the following inequality has been shown in [39], labeled Lemma 8 (p. 469) and (56) (p. 479): 
Lemma 2.4. For ([2.31) taken with h{x,t) = —{ut + Uux)ext, we have 


for 0 < s + r] < 1/2 and t > t*{U, H). 


(2.9) 


Remark 2.4. The above bound was critical in the previous analyses of this problem which made use of inertial 
terms —aAutt and strong damping —k 2 Aut [Hi or thermal smoothing [391140] . 


2.3 Nonlinear Plates with Delay 

In this section, the full flow-plate problem is reduced into a delayed plate problem. This delayed plate is 
subject three conflicting actions contributed by flow: 

1. The delay term itself: q'^(t). 

2. The term Uux in the plate equation—nonconservative, potentially leading to chaos. 

3. The term ut (with correct sign) in the plate equation, generated by flow. 

Above, (3.) is stabilizing. But (1.) and (2.) are “bad” and do not stabilize the plate. It is also known that 
a delay term can produce destabilizing effects on the dynamics, unless the damping (fc -|- l)ut is sufficiently 
large m- This is one contributing factor for our need of a large damping— k value. 

On the other hand, the term can typically be dealt with via a large static damping |32j . Here, we do 
not assume any static damping for the overall model . Our challenge, then, in this paper is to show that this 
potentially chaotic term is “thrown” to the equilibria set, and that the dynamics converge to it. 

The fundamental work in m focuses on the reduction of the model in (HI) to a delayed plate equation: 

Theorem 2.5. Let the hypotheses of Theorem M. 1\ he in force, and {uq, ui; (po, (pi) G ifg (H) x L 2 (Ll) x FI 1(R3 )x 
L 2 (Ri 4 _). Assume that there exists an po such that (po/x-) = (piii^) = 0 for |x| > pQ. Then the there exists a 
time L^{po, U,Ll) > 0 such that for all t > the weak solution u{t) in 11.1\) satisfies the following equation: 

utt + A^u + kut + f{u) =po- {dt + Udx)u- q'^{t) (2-10) 

with 

g“(t) = — [ ds [ d9[MgUext]{x — {U + sm6)s,y — scos6,t — s). (2-11) 

27r Jo Jo 

Owing to this reduction we can choose to study a general delayed plate model (as was done in [16]) 
to produce maximal generality with respect to our results for the delay system (which can be viewed as 
independent of our results for the full flow-plate dynamics): 
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( 2 . 12 ) 


{ Utt + + koUt + f{u) = po — q{u*) + Lu in il 

u = dvU = 0 on r 

u(0) = uo, ut(0) = Ml, M|(-f,o) = q 

The linear operator L : i^o(^) corresponds to non-dissipative lower-order terms : Lu = —Uu^ 

and = k + \. We have utilized the standard notation that m*(s) = u{t -I- s), so w* denotes u{s) for all 
s G {t — The operator q{u^) represents a general delay type operator (depending on u‘), of which (/“(t) 

is a specific example. (This notation is consistent with |16).i 

We will also need the following formula for the time derivative of (7“(t) appearing above in (I2.11|) : 


^27r 1 

dt[q^m= —[M^uUti^{U,9,0),t)d9 


(2.13) 


27r 


[M^uUt{K{U,9,t*),t-t*)d9 


0 ^0 
t* p27T 


{U -I- sin6») — [MgUa]e^t (x([/, 6», s), t - s) d9ds 




The proofs of Theorems 13.11 and 13.21 cited below are based on an analysis of the delay evolution (rt,H), 
with H = i?o(n) X L 2 {LI) X L 2 (—t*, 0; iJg (fl)) , corresponding to the reduced plate with given data ccq € H 
we have that Tt{xo) = {u{t),ut{t);u*) with Xq = {uo,ui,r]). The norm is taken to be 


||(m,u; 77)11 ^ = ||Am ||2 -p ||u|p + J \\Ar]{t + s)\\'^ds. 


In application, we will consider an initial datum yo G Y corresponding to the dynamics in (II.1[) (the full flow- 
plate dynamics) at time t = 0. We employ the reduction result Theorem [231 and we may consider the “initial 
time” {t = to) for the delay dynamics corresponding to any time after the reduction time L^{po, U, LI) above. 
At such a time, the data which is fed into ()2.10p is xq = {u{to),ut{to),u^°), where this data is determined by 
the full dynamics of (11.11) on {to — t*,to). (That (11.11) is well-posed on Y and (|2.12l) is well-posed on H (see 
[16]) allows us to move back and forth between the two descriptions of the system.) 


2.4 Attractors 

A global attractor A for a dynamical system (TjjH) is an invariant, compact set which is uniformly attracting, 
in the sense that if S C H is a bounded set, we have that 

lim d-ii{Tt{B),A) = 0, 

t—¥00 

(where we have utilized the Hausdorff semi-distance in H). 


3 Supporting Results 

For self-containedness of exposition we now provide a summary of the key results of m which will be used 
critically in this treatment. 

3.1 Plate Dynamics 

The first theorem we present is the main result of [16], given there as Theorem 3.4. 
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Theorem 3.1. Consider any nonlinearity of physical type. Suppose 0 < U ^ 1 and po G L 2 {^). Take any 
k>0 Then there exists a compact setU C Hgift) x L 2 {pfj of finite fractal dimension such that 

lim dy ({u{t),ut{t)),u) = lim inf (||M(f) - vo\\l + \\ut{t) - vif) = 0 

t—)-00 \ / t—>00 

for any weak solution (u,ut]4),(ft) to fi.il) . with initial data (mq, mi; ^ 0i </>i) G Ypg, i.e., which are localized in 
—d'o(x) = '/'i(x) = 0 for |x| > Po for some po > 0. We have the extra regularity U C fl x 

and any plate trajectory (u,Wt) on the plate attracting set (namely with initial (uo,Ui) GlA) has the 
additional property that 

{u,uuutt) GC{R; x H^{n) x ^2(0)) . 

In the setting of [32] and the present analysis we are interested in obtaining a sufficiently large damping 
so as to guarantee strong convergence of full flow-plate trajectories. It is imperative here that there exists 
a uniform absorbing ball for the plate dynamics which is independent of the damping parameter (when it is 
sufficiently large). The following result is given Lemma 4.6 in [32] (there applying to (11.131) for (3 > 0). 

Lemma 3.2. With the same hypotheses as Theorem \S.l\ above, there is an absorbing set B C Hq[VI) x L 2(^) 
which is (i) uniform with respect to the initial data yo (though it may depend on po ) and (ii) uniform in the 
damping parameter fc > fc* > 0, for some fc*. 

Remark 3.1. When we include the effects of static damping—of the form fiu, as in (11.131) — in the plate 
dynamics, we also have the existence of the attractor (and corresponding absorbing set) whose size also can 
be made independent of /3 > 0. 

The results in [T3] perform analysis on the abstract delay equation f2.12l) on H, including energy methods 
which show the existence of a true global attractor (as well as an absorbing set D '^) for the the delay 
evolution (rt,H). These results depend on taking ko > 0 (where ko is the damping coefficient in (12.121) '). 
(Note that, for the specific delayed plate of interest here (12.101) . we have ko = k + 1, and thus we could take 
k G (—1,00].) Then, to achieve the results in Theorems 13.11 and 13.21 projections are taken which suppress the 
delay component of the evolution. 

Remark 3.2. In what follows, the notation ^ GL correspond to the global attractor and (a fixed) uniform 
absorbing set for the dynamics (Tt, H). The notation U C B will refer to the attracting set and absorbing ball 
for the projected dynamics (u,ut) of (St,V) (delay and flow components dropped). 

3.2 Weak Stability 

Given the existence of a compact global attracting set for the plate dynamics, as well as the finiteness of the 
dissipation integral fCorollarv ll.6L we have the following theorem (given as Theorem 7.1 [32] and applying 
for (11.131) for /3 > 0): 

Theorem 3.3. Assume k > 0. For any initial data yo = {uo,ui] (fo, (fi) G Y and any sequence of times 
—>■ 00 there is a subsequence tn- and a point y = (u, w; (j), if) so that: 

St„. (yo) V- 

2 . ||M(tnj ) - wlh.n-t 0. 

3. ||Mt(t)||o —t 0, t —>■ 00, and hence w = 0. 

j. Along the sequence of times tn^ 00 we have 

sup ||M(tnj + t) - u||2-5,n-t 0 for any fixed c> 0, S G {0, 2). (3.1) 

rG[-c,c] 


In view of the above result, it is clear that the general, novel mathematical challenge in resolving the flutter 
problem is to improve the convergence of the flow component of the dynamics from weak to strong. This is 
a highly non-trivial task, taking into consideration the lack of compactness properties of the dynamics, along 
with the lack of dissipation acting on the flow equation. 
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4 Prom Weak to Strong 

We present (in a reorganized way) the key past results in [321 SO] , which yield supporting results that can then 
be applied in the proof of our main results here. Additionally, we restate in the present context certain ideas 
and results which were critically used in [32] . 

4.1 Hadamard Continuity 

The key results below address the continuous dependence property of the dynamics, namely, that St{-) is 
Hadamard continuous on Y (and Yg on finite time intervals). This was first addressed within the proof of 
Theorem 3.1 [32J pp. 1963-1964], and the proof of both Lemma I3T] and Theorem l4.2l can be found in [33]. (In 
fact, the Hadamard continuity results below are analogously true for von Karman dynamics.) 

Lemma 4.1. Let 0 < U < 1 and assume po € L 2 {Ll). Assume yo = {uo,ui](j)o,4>i) G Y. Consider the 
dynamics generated by generalized solutions to (HH (or (11.1311 if (3 > D), denoted St{yo)■ For any k > 0, /? > 0, 
and for any fixed T > 0, we have 

St{y^) ^ y{t) = Stiyo) G Y, 

when y™ —>■ yo G Y, uniformly for all t < T; this convergence may depend — possibly — on T > 0. 

Remark 4.1. The corresponding delay system (Tj, H) exhibits the Hadamard continuity property on any [to, T] 
as well. This follows from Lemma I4d1 and Theorem l2.5l Moreover, in m the continuous dependence property 
is shown directly from energy estimates on the delay system (12.121) . 

Remark 4.2. At present, it is not clear that the full dynamics for (EH are global-in-time Hadamard continuous. 
Whether this is the case is related to the decay rate of ||Mt(t)||L2(n) Ol indeed, to obtain this global-in-time 
property it is sufficient for ut G Li(r, oo; L2(H)). With fc > 0, we do have ut G ^2(0, oo; L2(H)); however, 
improving L2(0,oo) to Li(0;c») in general requires very substantial work. See Remark 14.101 

If we consider the addition of large static damping (see Remark 13.Ill , as well as large viscous damping, we 
see that the Hadamard continuity can be improved to global-in-time (this was one of the primary results in 

[32]): 

Theorem 4.2. Let 0 < U < 1 and assume po G L 2 {Ll). Assume yo = (mq, ui; ((q, (j)i) G Y. 

Consider the dynamics generated by generalized solutions to (jl.l3l) (i.e., with P > 0), denoted St(yo). 
Assuming that both damping parameters are sufficiently large k>kc>0, P> Pc >0, the semigroup St{-) is 
uniform-in-time Hadamard continuous, i.e., for any sequence y(f —> yo inY and any e > 0 there is an M so 
that for m > M 

sup||S't(y™) - S't(yo)||Fp < £• 
t>o 

This theorem and its proof are given in Section 7 of [32] . Theorem 14. 2 1 will be used below in an auxiliary 
fashion on a decomposed portion of the full flow-plate dynamics which resembles (11.1311 (i.e., has static damping 
added)—Theorem 15.71 

Remark 4.3. The minimal damping coefficients kc and Pc for the result above depend on the invariant set 
for the delayed plate dynamics, which itself depends on the loading po a-nd b, as well as the domain LI and 
other physical constants such as U, po, but is independent on the particular initial data of the system. 

4.2 Bounds in Higher Norms 

We now state and prove the critical supporting theorem for this treatment (and also [32]) Theorem 14.31 

For finite energy solutions to the flow-plate system, a uniform-in-time bound on plate solutions 
in a higher topology yields the desired convergence to equilibria result (as in Theorem 11.81) in the 
topology Yp for any p > 0. 

In other words, beginning with finite-energy solutions, if we further know that the plate dynamics are suffi¬ 
ciently regular, then compactness properties can be transferred to the flow dynamics. This result is independent 
of the particular structure of the plate dynamics—only global-in-time bounds in higher norms of the plate 
solution {u,ut) are needed; indeed, Theorem 14.31 follows from the structure of the (linear) coupling via the 
Neumann flow condition, written in terms of [w* -|- Uux\eyA^ and appearing in Theorem 12.31 
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Theorem 4.3. Suppose 0 < f/ 7 ^ 1 and po € L 2 {^), and take any fc > 0. Let (u, (j)) be a weak solution to (11.11) 
(or (11.131) if P > 0) with finite energy (flow-localized) initial data in yo £ Yp^. If there is a time T* so that 

sup {\\u{t)\\l + \\Aut{t)\\l + \\utt{t)\\l} <Ci, (4.1) 

[T* ,00) 

then for any sequence of times —>■ +c», there is a subsequence of times and a point y = (u, 0; cj), 0) with 

(u, 4>) € Af so that 

= 0 

for any p > 0. This implies that the result of Theorem \1.8\ then holds. 

We present this theorem as an independent result, in the spirit of what is used in [321 EO]- However, a 
word of caution: the assumed bound SH) will be valid only for “certain” solutions to the plate problem, and 
certainly not for every weak solution. This is in striking contrast with |391140] where the smoothing properties 
of thermoelasticity provide the additional boundedness of all plate solutions in higher topologies. Here, we can 
obtain the global-in-time bounds in higher topologies (as in (14.11) below) by considering smooth initial data 
for the flow-plate system and propagating this regularity—see Section 14.31 and 14.41 below. Alternatively, a key 
point of our proof below relies on a decomposition of the plate dynamics (see Section 15.21) in which one of the 
decomposed portions of the plate dynamics will satisfy a bound like that above in (|4.1I) (and the other will 
be exponentially stable). The non-decomposed plate dynamics will not necessarily have such a global-in-time 
bound. See Remark 14.41 below for further discussion of how Theorem 14.31 and its proof are utilized later. 

Proof of Theorem We first note that the solution to (12.21) (via classical scattering theory—see [321110] and 
references therein) is stable, so 

^ ( 0 , 0 ) 

in the local flow energy topology Yfi. 

We now utilize the assumption (SID, along with the bound in Lemma l2.4l to bound, uniformly, the flow in 
higher norms. Specifically, we note that by assumption 

u G C^{[T\ 00); L2(H)) n C\[T*, 00); H^iTl)) n C([T*, 00); n ( 4 . 2 ) 

(Below we utilize the notation 7Y®(iJ''(H)), meaning that regularity refers to the time variable. ) 

This implies, in particular, that u G iJ^(L 2 (H)) D D L 2 {H^iLl)). With u G L 2 {H^{Tt)) D 

H^{H‘^(Vt)), via interpolation [51133]. we have 

u G ai G [0,1]. (4.3) 

Similarly, from ut G L 2 {H^(fil)) D {L 2 (fil)) we obtain 

ut G a2G[0,l] (4.4) 

To apply inequality in Lemma 12.41 with s = 0, ,5 > 0 we will take ai = /3 = 02 


4 — 2ai > 2 -I- ai; so, 2 > 3ai and 2 — 2 q;2 > I + CI 2; so, 1 > 3a2. 


Hence, for ry < 1/3 we have: 


where Ci corresponds to the global-in-time bound on plate solutions in the assumption of the theorem. Then, 
applying the compactness criterion in Lemma 12.21 we have shown: 

Lemma 4.4. For any sequence of times —>• -foo, there is a subsequence (also denoted tn) and a limit point 

(u, 0; 4),ip) €Y such that: 

{u{tn),Ut{tn)](){tn)** 

in Y as n ^ 00. (Note: we know, a priori, that ut ^ 0 in L 2 {Ll)—see Theorem \3.tA ) 
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Now, to further characterize the flow limit point, we return to the formula in Theorem 12.31 

Lemma 4.5. The limit point tp gY in Lemma \4.4\ above is identified as 0 in L 2 {Kp) for any p > 0. 

Proof of Lemma EH The flow solntion {(j)*{t), (flit)) tends to zero (again, we emphasize that this convergence 
is in the local flow energy sense); hence, we have that (j){tn) </>■ To identify if with zero, we note that (by 

assumption) ut{f) G This allows to differentiate the flow formula in Theorem l2.3l in time (as was done 

in m and [30]) to arrive at that in (12.71) . 

For a fixed p > 0, we multiply <j)l*{x,t) in (12.71) by by a smooth function ( € C^{Kp) and integrate by parts 
in space—in (lO) we move dx onto C in the third term and the partials dx, dy from Mg onto C in term four. 
This results in the bound 


Wt\C)L2{K,)\<C{p) sup |lMt(t-T)||o.n||Clli.iCp (4-5) 

From this point, we utilize the fact that Ut{t) — >■ 0 in L 2 {Lt) when t —>■ oo (Theorem 13.3|) . and hence 
((^“(t), C)ifp —)> 0 and — >■ 0 in L 2 {Kp) for any fixed p > 0. This gives that (pt{t) ^ 0 in L 2 {Ll), and we 
identify the strong limit point ^ with 0 in L 2 {TL) so that {u{tn),ut{tn)]<p{tn)**,(pt*{tn)) — >■ (w, 0;((i, 0) in Y as 
n oo. . □ 


Now we show that the limit points obtained above satisfy the static problem in a weak sense. 

Lemma 4.6. The pair {u,(p) as in Lemma \4.4\ and Lemma \4.5\ satisfies the stationary problem (|1.22[) in the 
variational sense. 

Proof of Lemma \4.6\ We begin by multiplying the system (HD) by p G C'(j“(f2) (plate equation) and ip G 
((("“(Ri).) (flow equation) and integrate over the respective domains. This yields 


{utt, p) + (Au, Ap) + fc(ut, p) + (/(u), p) = (po, p) - {tr[(pt + Ucpxfp) 

{(ptt,tp) + U{(ptx,tp) + U{(pxt,tp) - = -{V(p,\7tp) + {ut + Uux,tr[tp]). 


Now, we consider the above relations evaluated at the points (identified as a subsequence for which the 
various convergences above hold), and integrating in the time variable from to tn + c. Limit passage on 
the linear terms is clear, owing to the main convergence properties for the plate component in the Theorem 
1331 The locally Lipschitz nature the nonlinearities of physical type allows us to pass with the limit on the 
nonlinear term (this is by now standard, [13]) ■ We then arrive at the following static relations: 

f {Au,Aw) + {f{u),w) = {po,w) + U{tr[^],Wx) 

\ (V<P, V^) - = U{ux, tr[iP]). ^ ’ 


This implies that our limiting point (it, 0; cp, 0) of the sequence 

{u{tn),ut{tn);(pitn),(pt(tn)) is in fact a solution (in the weak/variational sense) to the static equations, i.e., 
it is a stationary solution to the flow-plate system HU). We have thus shown that any trajectory contains a 
sequence of times {tn} such that, along these times, we observe convergence to a solution of the stationary 
problem. □ 

The above lemmata complete the proof of Theorem 14.31 

□ 


Remark 4.4. In the proof above we have shown how a bound on the plate dynamics in higher topologies results 
in the desired convergence to equilibrium result for the full dynamics. We pause to point out that this will 
be used in two different ways: below we will show that with smooth initial data for HU, via a propagation 
of regularity result, such a bound can be obtained. Alternatively, in the proof of our main result (in Section 
15.5|) . we will use the proof of Theorem 14.31 for decomposed systems (15.251) and (15.271) . We note that the result 
of Theorem 14.31 is not strictly dependent on the structure of the plate equation, but rather the linear nature 
of the coupling via the Neumann condition for the flow (12.31) . Hence, if we add change slightly the structure 
of the plate equation, the result of Theorem 14.31 will be unchanged, so long as (i) the modified plate dynamics 
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still converge in the sense of Theorem l3.3l and (ii) the requisite bound in Theorem l4.3l is obtained. Specifically, 
in what follows we will consider a general system: 

{ utt + (k + l)ut + = pq — UdxU — + F{t, u) in Q. 

u = d^u = 0 on r (4.8) 

u(to) = UO, Ut{tQ) = Ui, = ?7, 

where the data will be “smooth” (as in the proof of Theorem l5.6l) . or possibly zero (as in the proof of Theorem 
I5.10|) . In both of these cases the nonlinearity F(t,u) will be broader than / = /(u), depending on a given 
function^ owing to the decomposition of the dynamics considered below in (15.1|) and (15.2|) . 

4.3 Propagation of Regularity 

In order to utilize the key technical supporting result in Theorem l4.31 we must show a propagation of regularity 
result for the dynamics. Specifically, we must show that for “smooth” initial data, we have “smooth” time 
dynamics; this is Theorem 14.71 Such a result was correctly stated in the [3H (2014)], however the proof was 
unclear. This has been updated in a subsequent version [3H (2015)] of the same reference. We include a proof 
here in the case of Berger dynamics (with the clear proof for von Karman dynamics in [321 (2015)]). 

Remark 4.5. To show infinite-time propagation of the plate dynamics, we rely on the full flow-plate dynamics 
to achieve propagation of the initial regularity on any interval [0,r*]. Once this is achieved, we may work 
on the reduced delay plate (after sufficient time has passed) and utilize sharp bounds to obtain regularity of 
plate trajectories on the infinite horizon (propagation on [T*,oo)). This part of the argument depends on (i) 
uniform exponential decay for the nonlinear plate equation with large static and viscous damping [30] . and 
(ii) specific properties of nonlinearity; (in the von Karman case one applies the sharp regularity of Airy stress 
function [T?] p. 44]. 

Theorem 4.7. Consider the dynamics {St,Y) corresponding to (II.ip . Consider initial data ym € ^{T) such 
that y™- G BrIY) and take k > 0. Then we have that for the trajectory St{y'^) = (w'"(t), M™(t); (/)'"(t), (^]"(t)) 

G Cl (0,r;ij2(O) X H\Rl)) , 

for any T, along with the bound 

sup {WAuY^Wl + \\u^\\l} < Cm,T < oo. (4.9) 

te[o,T) 

Additionally, if we assume the flow initial data are localized and consider the delayed plate trajectory (via the 
reduction result Theorem ] 2. 5\} . we will have: 

sup {\\AuT\\l + \\nTAn}<Ci<oc. (4.10) 

tG[0,oo) 

Then, by the boundedness in time of each of the terms in the (12.101) , we have 

sup \\A'^u^{t)\\o < C 2 , 

t^[0,oc) 

where we critically used the previous bound in (14.101) . In particular, this implies that, taking into account the 
clamped boundary conditions: 

sup l|u'"(t)l|4,a < C3. (4.11) 

tG[0,cso) 

Each of the terms Ci above depends on the intrinsic parameters in the problem (including the respective loading, 
which depends on the nonlinearity being considered). 

Remark 4.6. The superscript m in the above theorem serves to indicate that “smooth” data (in the sense 
above) will be used to approximate given finite energy data. 

Remark 4.7. The parameter ,5 > 0 is inert with respect to this result. That is, Theorem 14.71 holds for (11.11) 
(/3 = 0) and (11.131) (/3 > 0). 
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For this proof we consider no static damping (/3 = 0 in (11.131) ') and any k > 0 fixed. Additionally, this 
propagation result holds for both the Berger and von Karman nonlinearities [32] . 

The proof of this Theorem 14.71 is based on two steps labeled STEP 1 and STEP 2. We first prove finite 
time propagation of regularity (STEP 1). Following this, we shall propagate regularity for all times uniformly 
(STEP 2) without making assumptions on the damping coefficients. This is motivated by the fact that for 
the infinite time propagation we must use the delay representation for plate solutions, which is valid only for 
sufficiently large times. Thus, the regular initial condition required to do so is obtained via the finite time 
propagation of regularity for the full flow-plate model. 

Proof of Theorem 

4.3.1 STEP 1 

To prove (14.101) above (which then implies (14.111) 1 we will consider the time differentiated version of the entire 
flow-plate dynamics in dn]), which is permissible in Stiff). We label u = u™ and 4> = (/)™, multiply the plate 
equation by Ut and the flow equation by (/)(, and integrate in time. We consider the ‘strong’ energies 

E;i{t) = \[\\^u[t)r + \\Mm\ 

Ejiit) =l[\\mt)\\ii - 
EUt) = 2U{Mt),^{t))n 

£^t) = Eli{t) + + Et^,{t), E^(t) = Eliit) + E^fi{t) (4.12) 

From the standard analysis (done at the energy level) we have that [Ill [44] 

\Et^M < ^llVCWf + j\\u7M\l < <5||VCWf + ellAurWf + c.,s\WTit)r (4.13) 

Using the boundedness of trajectories at the energy level sup^^o ll'*^r(Ollo < Cr, we then have that for all 
e > 0 

\EUt)\ < + K,, (4.14) 

and this results in the existence of constants c,C,K >0 such that 

c-E%t)-K<£%t)<Ci:,'^{t)+K, Vt>0. (4.15) 

Then the energy balance for the time-differentiated equations is: 

£^{t) + k [ \\ut{r)fdT=£^{0) + B{t), (4.16) 

where B{t) represents the contribution of the nonlinearity in the time-differentiated dynamics: 

B{t)^ (4.17) 

Using (14.151) . we have 

W{t) <CW{0) + K+\B{t)\ (4.18) 

We note immediately that if we can bound B{t), we will obtain the estimate in (14.101) . Continuing, we use the 
key decomposition for f{u) = [b — ||Vm|P]Am: 

(1 = 2(Vu-,<)(Au™, 0 + [||Vu-|p - b]iAuT,ufi) (4.19) 

= j^Qiit) + Pi{t) (4.20) 
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where 


(4.21) 

(4.22) 


Qi{t) = - - b]\\VuT\? - {Au^,uTf, 

Pi{t)= -3(AM™,0||Vit™||2. 


Thus we have 

\B{t)\ < Clll [||Vu-|p - b] ||V<IP - {Au^,uTno\\uT\\o\l+ j\{^u^,uT)\ llVuriP^r} 

< e\\uTit)\\j + C.imuTW^ + C(E«( 0 )) + Cr f WuTWjdr 

Jo 

We again appeal to the hyperbolic estimate: 


(4.23) 


\uTmi<\\< 



\UttiT)\\ldT, 


which yields: 


\B{t)\< e\\uT{t)\\l + aA^%0)) + C,,R / [\\uT\\l + \\uTtro]d7 


< Wul 


‘(t)||2+C,.fl(E*(0))+C,.fl [ E*(r)d7 

^0 


(4.24) 

(4.25) 


Hence the above inequality, when used in conjunction with (14.1511 . (14.1611 . and Gronwall’s inequality, provide 
uniform boundedness of the higher norm energy E®(t) on any finite time interval [0,T], for any finite T. 


4.3.2 STEP 2 

Remark 4.8. Using the variation of parameters formula, and global-in-time bounds of delayed plate trajectories 
in the energy topology, one can obtain the bounds in (14.1011 and (14.1111 straightforwardly, however this requires 
large viscous damping (see (15.1411 1. Here we show below that these bounds (and the final result for strong 
stability for smooth initial data) can be obtained with minimal viscous damping and NO static damping. 

In order to obtain uniform-in-time boundedness for plate solutions (as in 14.101 and 14.1 Ill , more work is 
needed. The needed argument depends on the reduction in Theorem l2.51 where the flow dynamics are reduced 
to a non-dissipative and a delayed term in the plate dynamics. This reduction is valid for sufSciently large 
times—thus the result of STEP 1 applies to all the times before the delayed model is valid. 

We return to the delay representation of the plate solutions (Theorem I2.5|l and utilize our finite time 
horizon regularity from STEP 1 to obtain regular plate data at time T. Considering (u,ut) = we 

write 

utt + A^u -I- (fc -I- l)ut = ^ {[||Vm’"|P - 6]Au™} - UdxU- ^ {g(u'")} (4.26) 

which is valid for any t > max{T, t^}. 

We utilize an approach taken in m for the non-delayed plate. We will define a Lyapunov-type function 
for the time-differentiated, smooth dynamics so that we may use Gronwall, along with the finiteness of the 
dissipation integral (Corollary [T^ : 



lulrdr < -foo. 


Let us use some more convenient notation here: E(t) 



-t-||Mt|p]. Then define 


Q{t) = E{t) - Qi{t) + z/||u| 


where Qi and Pi are as above. 
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Let z/, /i > 0. We can choose z^(i?) large enough so that Q{t) is positive [T3]; we make this choice and then 
keep z/ fixed. We will now construct a Lyapunov-like function for the dynamics T^{xq) = u*) £ H 

with xq £ H. There are constants ao(z^) and ai(z^) so that: 


ao 


IIAMr + lluzir <QW<ai \\Mr + \\ut 


The parametric Lyapunov function is then: 

W{t) = Q{t) + e{ut,u) + [ [ E{T)dTds, e£(0,eo)- 

Jo J t — s 

The parameter /i will be addressed below. The value of Eq is chosen so that 

1 /■* 

-aoE{t) < W{t) < 2aiE{t) + fit* / E{t), V e £ (0, Eq), 

^ t — t* 

fi f f E{T)dTds < fit* f E{t). 

Jo Jt — s J t — t* 


and we note that 


We compute W'{t): 


Make the replacement: 


^W{t) = {Au,Aut) + {ut,utt) - + v{u,ut) 

+ e{uu,u) + e\\ut\\'^ + fit*E{t) - fi f E{t) 

Jt-t* 


utt = -A^u -{k + l)ut - —q'""' - -T:[f{u^)] - Uux- 

at at 


Then: 


—W(t) = {fit* -I- E - A: - l)||ut||^ -I- (z^ - e/c - l)(u, Ut) - ([g“ ]', ut) 
-{[f{u-^)y,Ut)-U{u,,Ut)-j^Qi 

{fit* - e)\\Au\\'^+e{[q'^'"]',u) - £{[f{u)]',u) - £U{ux,u) - fi [ E{t) 

Jt-t* 


(4.27) 


(4.28) 

(4.29) 


Thus: 

d f* 

-W{t) + {k + l)\\utf+e\\Auf+fi E{T)dT (4.30) 

at Jt-t* 

= {fit*+e)\\utf+fit*\\Au\f + P,{t) - E([/ («“)]',«) 

-e([ 9 “ ]',u)-([g" ]',ut) + {ii - e{k + l)){u,ut) - eU{u::o,u)-U{ux,ut), 

Via (14.271) . choosing fi and e sufficiently small, there is a 7(01, z/, k, fi, t*) > 0 so that for any e £ (0, Eq): 

^W{t)+j[E{t) + £’(T)dT] < c|Pi(t) - e([/('u)]',u) - E([g“™]',M) (4.31) 

- {[q'^"']\ut) + {ti-ek- l)(u,?It) - eU{u^,u) - U{ux,ut)'^ 

Lemma 4.8. There exists a 7 > 0 (with the same dependencies 7 above), and constant M (depending on R 
and the same quantities as j) so that: 

W'{t)+^W{t) < M[l + \\u\\'^E{t)] (4.32) 
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Proof of Lemma \4.8[ We recall (12.131) . 
energy level (u™, u™), 


So, along with the global-in-time priori bounds for the lower energy 


sup 

tG[to,oo) 




< OO, 


and 


= _ 2(Vu’", VOAit"" -h [6- ||Vu’"inAu™ 
= 2(Azi'",M)Au'" + [b- ||Vu'"|p]Au, 


we obtain: 


|Pi|<C||Au-|| INI llVuN (4.33) 

\{[f{u-^)r,u)\<C{R)[\\u\f + \\u\\ IIANI] (4.34) 

(v - ek - l)|(u,?It)| < (5||ut|P + Cs{v,£,k,R) (4.35) 

£U\{ux,u) + U\{ua:,ut)\ < ^[||Mt|P + IIAmN] -h Cs{£,U,R) (4.36) 

e|([<7“"‘r,«)|<C(e)(^|K(t)||2 + ||u™(t-t*)||2 + + (4.37) 

-l-(5||u||2 -l-C5(e)||M||^ 

<5E{t)+Cs{£,t*,R) (4.38) 

And, noting the explicit estimates for the term q'f performed in [121 (6-1)], we also have 

{dt[q-^ln,) < <5||zZ*|p + Q(||u-(t)||i,o + \\u^{t - t*)\\ln (4-39) 


Then we note: 


Nli<<5||u||l + C,(A-)||u||BC. + e||u||i 


where ||m ||4 is to be calculated from an elliptic equation in terms of || (u, Ut)\\Ypi- Indeed, we consider biharmonic 
problem with the clamped boundary conditions. 


with 


This gives 


We have 


= -ut -{k + l)u+ [||Vu™|p - b]Au^ - Uuff - 

u™ = Nw™ = 0 on dkl. 

IN^IU < C\\ut -{k + l)u+ [||Vm”^|P - b]Au^ -u^- 
<Cn|(u,zlA||v, + ||u’"||i + ||u™||3 + ||u™||21. 


\\un\i<cmEit) + M- 


([g“'"r,iit)<l|i7t||o 




< S 


Eit) + / E{T)dT 


Cs{t*,R) 


(4.40) 
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Combining these estimates, we return to (j4.31l) : 

^W{t)+^[E{t)+ ^ E{T)dT] = C'|Pi(t) -£([/(«’")]', m) -e([g“’"]',M) - m*) 

+ {v - e{k + l)){u,ut) - eU{uj:,u) - U(ux,ut)'^ 

< 6[E(t) + f E{T)dT] + Cs{e,v, R)\\u\f E{t) 

Jt-f- 

+ Cs{R, iy,e,t*,U,k) 

Scaling S appropriately and absorbing terms, we have: 


d_ 

dt 


W{t)+j[E{t)+ [ E{T)dT]<M + C\\u\\‘^E{t). 
Jt-t* 


Recalling (I4.27L we have some 7 > 0 such that: 


—W{t)+W{t) <M + C\\ufE{t). 


Noting (I4.27|) . we have: 


dt 


W{t)+^W(t) < M + CiWuW^Wit), 


and utilizing the integrating factor we see that 

d 


— + Cie~^*\\u\\'^W{t). 


Integrating, we have: 


e-fTw{T) < e^*°W{to) + ^ 

1 L 


e7T _ g7to 


Cl / e^^\\u{t)\\'^W{t)dt, 


to 


and thus 


W{T) < ^ 

7 L 


< ^ + e-^(T-to) 

7 


Wito) - 


I _ ^-l{T-to) 

M 


■Cl / \\u{t)\\^W{t)dt 


'to 


7 


+Ci / \\um^w{t)dt. 


Let 


(T) = ^ + 


7 


'to 


M 

W{to) - 

7 J 


and note that a{T) > 0 for T > to. Then, using Grownall, we have: 


(4.41) 


(4.42) 


□ 


W{T)<a(t)+J^ a(t)||u(f)||^ exp ||u(s)||^ds| dt. 

Note that a{T) > 0 and a is bounded, and thus, owing to the finiteness of the dissipation integral 

^00 nOQ 

/ ||M|pdt= / ||M™||Qdt <+00, 

J to ''to 

we obtain global-in-time boundedness of W(t) . Then, using (I4.27P , we have global-in-time boundedness of 


E{t) = -[\\Au 


riP + iKiP] 
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Via elliptic theory (as utilized above), we also have ||M'"(t )||4 globally bounded in time. We note that tV(to) 
has the requisite regularity by STEP 1 above (where we have shown propagation of regularity on the finite 
time horizon for the full flow-plate dynamics). 

This concludes the proof of Theorem 14.71 □ 

Remark 4.9. As mentioned above, the proof steps can be easily adapted to the case of the von Karman 
nonlinearity using the sharp regularity of the Airy Stress function [141 p. 44]. 

4.4 Strong Convergence for Smooth Data 

From Theorem 14.71 we see that for the Berger dynamics (and, in fact for von Karman), in the presence of 
minimal viscous damping (only requiring fc > 0), we have a propagation of regularity result on the infinite¬ 
time horizon. From this bound, for initial data j/'" G .^(T), we can immediately apply Theorem 14.31 This 
yields one of the primary results of [35J Theorem 2.1]: 

Theorem 4.9. Consider pm S Stiff) fl Br{Y) C (it^ fl Hq){CI) x x i7^(]R^) x C Y. Suppose 

that the initial flow data are supported on a ball of radius po (as in Theorem \2.5\) . Suppose k > 0. Then for any 
sequence of times > oo there is a subsequence of times tn^. identified by tk and a point y = (u, 0; (j),Q) G Y 
, with (m, (j)) G Af, so that 

Il-S'tjy'") - yllvp-)■ 0, tk^oo. 

This implies the primary result: strong convergence of (u{t),(j){t)) to the equilibria set JV for initial data 
yoGS(T)nYp,. 

Remark 4.10. We note that the situation when stability to the set of equilibria is achieved for smooth initial 
data is typical in the analysis of strong stability for hyperbolic models with damping. The next step is to 
extend the desired property via density to all initial data residing in the finite energy space. This can be done 
(see, e.g., 0 ) via a uniform-in-time Lipschitz (or Hadamard) property of the semigroup. As we describe in the 
next section, this is precisely the key issue to contend with in the analyses in [32] (where large static damping 
was used in conjunction with large static damping) and here, where a fortuitous decomposition allows us to 
consider only large viscous damping. 

4.5 Strong Convergence for Finite Energy Data—Technical Difficulties 

With the technical background established in Sections [21 14.31 we now address the specific difficulties involved 
in showing a stabilization to equilibria result without assuming either (i) terms of the form kiAuu — k 2 Aut 
(as in [nmi]), or (ii) exploiting parabolic effects in a thermoelastic plate IMlIlDj- In both cases, as well as 
in m, the key task is to first show compact end behavior for the plate dynamics. This requires the use of a 
reduction of the full flow-plate dynamics to a delayed plate equation (Theorem [23] above), at which point one 
may work on this delayed plate system. In both case (i) and (ii) above the ultimate character of the nonlinear 
component of the model is compact —owing to the fact that parabolic smoothing and rotational inertia both 
provide Vut G L 2 {Cl). This permits compactness of the boundary-to-flow (Lopatinski or hyperbolic Neumann) 
mapping. The results in m were the first to show that dissipation could be harvested from the flow in order 
to show ultimate compactness of the plate dynamics without utilizing the compactness of this mapping. 

We note that in [40] the key to the stabilization result lies in a compactness criterion (given as Lemma 10 
in |39| p. 472]) for flow trajectories (given here as (12.41) 1. The flow bound in higher norms is obtained via the 
thermoelastic character of the structural dynamics. 

In [32] . the major contribution was the ability to circumvent the seeming lack of natural compactness in 
the dynamics (for von Karman or Berger) by utilizing large static and viscous damping. The requirement for 
large static and viscous damping was needed to show exponential decay of plate trajectories, which resulted 
in uniform-in-time Hadamard continuity of the full flow-plate dynamics (see Theorem 14.21 and Remark 14.101 
above). With this, and the desired result holding for regular initial data fTheorem I4.9[) . an approximation 
argument was used to pass convergence properties from regular data onto finite energy initial data. The 
principal result there ( [321 Theorem 2.3]) is the analog of Theorem 11.81 however, it requires the addition of 
large static damping. 
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Theorem 4.10. Consider (11.131) . Let 0 < [/ < 1, and assume po G L 2 (Ll). Then there are minimal damping 
coefficients kmin and Pmin so that for k > kmin > 0 and j3 > Pmin > 0 any generalized solution {u{t), (j)(t)) to 
the system with localized (in space) initial flow data in Yp has the property that 

+ ll“tWllL(o) + Mt) - ^\\m{K,) + \\Mt)\\l,iKp)} 

= 0, for any p > 0. 

As in Theorem O here, the minimal damping coefficients k^i^ and /3mm depend on an invariant set for 
the plate dynamics, which itself depends on the intrinsic quantities in the problem. 

Though the results of [32] stand as the culmination of the treatments [mm, and critically use the results 
of [TU [39l [40] , static damping was necessary to contend with lower order terms (due to the non-dissipative 
nature of the reduced plate dynamics (12.121) 1. but seemed physically anomalous. The key technical point in 
this treatment is the elimination of this need for large static damping imposed on the structure. We seek the 
physical result which says that for large viscous damping only (size dependent upon the intrinsic parameters 
of the problem) the full-flow plate dynamics stabilize to the stationary set. Our approach here is based on a 
decomposition of plate dynamics into an exponentially stable component and a smooth component. To show 
smoothness of this latter component in the decomposition, we require working in higher topologies. 

5 Proper Proof of Main Results 

In this section we prove Theorem 11.81 The proof follows through several steps and critically uses the results 
and arguments presented up to this point. 

5.1 Outline of the Proof 

For the reader’s convenience we now give a summary of each of the key steps in the proof: 

In Section 15.21 we consider a decomposition of the delayed plate dynamics (12.101) : we consider a z and 
w component, such that z + w = u. We will decompose the nonlinearity in an advantageous way (for time 
sufficiently large), as well as add static damping fdz to the z equation and subtract it from the w dynamics. 
This allows us to use damping of the form kzt + (dz in the z plate dynamics, which will result in uniform 
exponentially stability for m- For the w component, we will consider zero initial data (the initial data 
being left in the 2 dynamics), giving (15.21) . The w dynamics, with zero initial data, will be smooth, having 
appropriate global-in-time bounds. 

Section 15.31 and 15.41 provide the proofs associated to the aforementioned plate decomposition and are the 
major technical contribution of this work. Section 15.31 follows |32] in utilizing Lyapunov methods to show 
exponential stability of z component of the plate decomposition; Section (|5.4p demonstrates the smoothness of 
the w component of the decomposed plate dynamics, and critically uses the null initial data in differentiating 
the dynamics and obtaining global-in-time bounds. 

In Section [5751 we decompose the flow according the plate decompositions in the previous step. This allows 
us to reconstruct a full flow-plate dynamics corresponding to z, and a full flow-plate dynamics corresponding to 
w, coupled through a mixed, nonlinear plate term. The (z, (p^) dynamics (corresponding to the exponentially 
stable 2 plate dynamics) will exhibit global-in-time Hadamard continuity (depending critically on large static 
and viscous damping), and thus we can approximate any trajectory by smooth data (this follows the general 
structure of the argument in [32]). The (w, (/)’") dynamics (corresponding to the plate dynamics with zero initial 
data) are smooth, and thus the results and techniques for smooth initial data can be re-applied (Theorems [12] 
andiTTi). 

5.2 Decomposition of Delay Plate Dynamics 

We consider the reduced, delay plate equation in (12.101) (see [TO] [H [16] ). The more general PDF system with 
delay (|2.12p is well-posed on H = Hq{Q) x L 2 {LI) x L 2 (— t*, 0; iLg(fl)) . We write Tt{xo) = (u, ut;u*) with 
xq = {u{to),ut{to)\u^°) G H. There is a compact global attractor ^ C H. [16], and we denote by C H 
a fixed absorbing ball for the dynamics of diameter R, such that C We note that for k > k^, > 0 
the size of the absorbing ball can be taken uniform in k. In what follows, we can reduce the above delay 
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plate dynamics (T^,!!) to those on the set ^ by waiting for absorption; this is to say that for any initial 
data ||(u(to),Mt(to);'a*“)||H < P there is a time of absorption tp such that Ttp{B-ii{P)) C Thus, by 
waiting sufficiently long (i.e., choosing to > tp —valid, as we are considering strong stability as t —>■ +oo), we 
can reduce dynamics corresponding to any initial data to dynamics with initial data chosen from We will 
denote this dependence by R —the diameter of 

Owing to the fact that q' is linear, we will apply the following decomposition of the solution (note this is 
a different decomposition than typically used for the proof of exponential attraction i.e., [25)1. 

Take {u{to), Ut{to)\u*’’) € C H. We write u = z + iu, where z and w correspond to the systems: 

ztt + (k + l)zt + + A^z + P(z, u) = — Uzx in V, 

^ P{z,u) = [b-\\Vu\\'^]Az 
z = dpZ = 0 on r, 

^z(to) = w(to), Ztito) = Ut{to), z\(^to-P,to) = 


Wtt + {k + l)wt + A^w = po + q'^ — Uwx — P{w, u) + /3z in fl 
^ = [6-||Vm|P]Au; 

w = dpW = 0 on r, 

^w{to)=0, wtito) = 0, = 0. 

Remark 5.1. Since we are taking null initial data for the the w portion of the decomposition, we note that the 
presence of z (via u = z -\- w) in the term P{w, u) is the driver of the dynamics. 

Given u, the z/w system is well-posed with respect to finite energy topology "H. Indeed, in light of the 
stability estimates shown below for the z dynamics (with u given) the analysis in [16] applies; from there, the 
w dynamics are well-posed as a difference w = u — z. We state this as a lemma: 

Lemma 5.1. Consider the dynamics in (15.11) at time to, with initial data (w(to), wt(to); m*”) G ^ C H. Given 
u £ C{[to — t*, r]; IJq (n)) the dynamics corresponding to (15.11) are well-posed in the sense of \10f . Then, 
considering w = u — z, with u corresponding to the solution in (12.101) on H, the problem (15.2p is well posed on 
H with {w,wt;w*) as a solution. 

Let 

Epiz{t))^l[\\Az\\^ + \\z,\\^ + (i\\zn 

For the decomposed plate dynamics, we will show the following supporting lemmas: 

Lemma 5.2. There exists ke{R) > 0 and (3e{R) > 0 such that for k > kg and j3 > jdg, the quantity Ep{z{t)) 
decays exponentially to zero. In fact, Tf{xo) = (z,zt;z‘) decays uniformly exponentially to zero in H for 
xq £ Bii{R) (with rate depending on k and (3). 

Lemma 5.3. There exists a kQ{R) such that for any k > kq all j3 > 0 the evolution (w,wt',w*) on H 
corresponding to (15.21) has that {w,wt,wtt) £ C{[to,oo)-, H^), where H'' = x iJ^(O) x L 2 {fl). 

We will need to utilize the following estimates on q, coming from |16j . It is shown that they are valid for q' 
as given in (I2.10p —see also ()2.13l) above. These estimates indicate a mechanism for compensated compactness 
acquired by the flow and translated into delay part of the plate system. 


||g“(t)|| < GG / ||Au(r)||dr, ||g“(t)||_i < Gt* / ||u(T)||idr, 

Jt-t* Jt-p 

\{qt{t),fj)\ <C^\\u{t)\\i +\\u{t-t*)\\i+ f ||u(r)||2dr|||7/;||i 

^ t — t* 


\qt(t)\\ < Ct* l|Au(t)||-f||Au(t - t*) 


f IK 

Jt-t* 


slllads 


(5.3) 

(5.4) 

(5.5) 
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The last inequality follows from (12.131) , after noting that Mg is the second order differential operator controlled 
by 11 Am 11 (due to zero Dirichlet boundary conditions imposed on the structure). Additionally, we note the 
following elementary bound: 


G{u{T))dTds <t* G{u{T))dT, 


where G is any positive functional. 


5.3 Stability of fl5.ip — Proof of Lemma 15.21 

We consider the Lyapunov approach which was used in [Ml [16] (based upon [13]) to show the existence of 
an absorbing set (and used in [32] to show Hadamard continuity of the full flow-plate dynamics). Define the 
Lyapunov-type function: 


V{T,^{xo)) ^ Ep{z) - (g(z‘),z) + (z*,z) + -||z|p 


(5.6) 


+ \\^z\fds + J ||A2:(r)||^drds), 


where T^{xo) = (z(t), Zt(t);z*) for t > to, and /r is some positive number to be specified below. Using Young’s 
inequality and (15.31) . we have immediately that 


c{^l)Ep{z{t)) < U(T/(xo)) < G{k)Ep{z{t)) + 


{z{i 


^dr 


(5.7) 


The constant c(/i) does not depend on either damping parameter. We now continue with the Lyapunov calculus 
by computing the time derivative of V{T^{xo))- Since the dynamics has been shown sufficiently regular for 
regular data, we can consider such regular initial data and perform differential calculus on a priori smooth 
solutions. Density and passage with the limit on initial data will allow to draw the same conclusion for hnite 
energy initial data. 

^U(T/(xo)) = ^Ep{z(t)) + {ztt,z) + llzill^ -h k{zt,z) 

- - {qt,z) 


^ it + J ||Az(r)|p(ir(is| 


> t — S 
|2 


= -A:||^t||^ + (^U+^-l)||Az||^-/3||z| 

WAzW^ds 


— fi\\Az(t — t* 


(5.8) 


-I- {q,z) - {qt,z) - U{zx,z) - U{z,c,zt) 

- {^{z,u),Zt) - {T{z,u),z) 

We twice made use of relation ([531): Ztt + (dz+{k + \)zt -I- A^z + T{z, u) = q^ — Uzx- We estimate from above; 
hrst we note that 


||J-(^,m)|| =11 [6-||Vm||2]Az || < C(i?,6 )||z||2. 


(5.9) 
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Then, using (15.31) . Young’s inequality, and compactness, we have: 

\{T{z,u),Zt)\ + \{H^.n),z)\< e\\^z\\^ + C,{RM[\\ztf + |k|P] 


\{q\t),z)\<e f ||Az|p + C,(r) 
Jt-t* 




\\Azit)r + \\Azit-t*)r+ / ||Az(r)||dr 

Jt-t* 

2 


+ Ce{t*)\\z„ 

U\{z,,z)\ + U\(z,,zt)\ < e[||A0||2 + ||^*||2] + C,{U)\\z\\^ 


We may choose sufficiently small, then e small, then fj (depending on R and fi) sufficiently large, and 
finally k sufficiently large (depending on /3 and R). 

This yields: 

Lemma 5.4. There exists ke and j3e so that for k > kg > 0 and j3 > Pe > 0.' 

jV(Tf{x)) < -c(fc,/3) + ||Az||2 + ||z|r + \\Az{t - t*)\f + \\Azir)\fdr^ , (5.10) 

where c{k,l3) —>■ oo as min{fc,/3} —>• oo. 

From this lemma, and the upper bound in (15.71) . we have for some 7(fc,/3) > 0: 

Ep{z{t)) < CV{Tf{xo)) < CV{xo)e~^\ t > to. 

This concludes the proof of the uniform exponential stability of the 2 portion of the decomposed dynamics on 
SS. 


5.4 Regularity of fl5.2p —Proof of Lemma 15.31 

Again, recall that we have chosen to large enough so that (u(to), Mt(io); Let T™(0) = {w{t),wt{t);w*) 

correspond to the (15.21) dynamics. First, we note: 

Lemma 5.5. The dynamics Tf^(0) (such that T(^(0) +Tf (xq) = Tt{xo), andu = w + z) are uniformly bounded 
on H. This is to say that: 

\\{wit),wt{t);w*)\\ii < C{R), V t > to- 

Proof of Lemma \5.5\ This follows immediately from the existence of a uniform absorbing set for the dy¬ 
namics Tt[xo) = (u(t), itt(t); it‘) on H, and the uniform exponential stability of the Tf{xo) dynamics. □ 

To continue with the proof of Lemma (lOl) . recall iF{w,u) = [fo — ||Vm|P]Aiu. Consider the time- 
differentiated w dynamics; let w = Wt- 

{ wtt + {k + l)wt + A^w = ^ {q{w*)} - {iF{w, u)} - Uwx + jdzt in ft 
W = a.zF=0 onF (5-11) 

Uj(to) = 0, Wt(to) = 0, ,ta) = 0- 

Recall, the term ||zt|| decays exponentially, and is thus bounded. We will now obtain a global in time 
bound on | |lc| 1 2 -b | \wt \ \ in what follows. We now invoke the well-known exponential decay of the linear, damped 
(static and viscous) plate equation for 

Wtt + {k + T)wt + Kw -f A^w = G. 
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(To note this decay one can, for instance, see the argument just given in Section 15.31 leading up to Lemma 
El) We consider: 


wtt + (A: + l)wt + Kw + - — {J'(w,-«)} - Uwx + fizt + Kw 


dt 


dt 


(5.12) 


= G{u, ZjWjWjW^). (5.13) 

Let w(t) = {W{t),wt{t)). For to (again, sufficiently large) we utilize the variation of parameters formula: 


f{t) 


Ypi 


< C 


'{to)) 


+ / 

Jto 




ikriioffi 


-Tiw,u) 


(5.14) 


+ 17||n;||i + K\\w\\o + /3|kt||o 


ds. 


Noting that uj —>■ oo as min[fc,iF] —>■ oo (see the analysis (I5.10|) l. The critical terms involve qf' and 
—J-{w,u), as the other terms are bounded by Lemma [S751 

For the term dt [<?’"] we note the bounds in (15.51) above (and [T^] for more details): 


MiT)\\<Ct* ||An;(r)|| + ||A«;(T-r)|| + 


k(s)||3ds 


(5.15) 


At this point, elliptic theory can be applied: a bound on ||rc ||4 is to be calculated from an elliptic equation in 
terms of || (W, Wt)||. Indeed, we consider biharmonic problem with the clamped boundary conditions. 

A'^w = —Wtt — {k + l)wt — u) — Uwj; — q'^ + fiz; w = d^w = 0 on dV,. 

This gives, via elliptic estimates (after accounting for the a priori bound in Lemma 15.51) : 

||w||4 < C\\wtt - ik + l)wt - T{w, u) - Wx - q"^ + l3z\\ 

< + (fc + Ijllwill + (1 + ||■u||?)||^c ||2 + C/llialli + (1 + t*) sup ||w(s )||2 + pc] 

< C[\\{wt,Wtt)\\Ypi + !]■ 

Thus we have: 


lk(A)||4 < C{R) [||(w,wO||yp, + l] , 


We also have: 

d 
dt 


J-{w, u) 


{Vu,Vut)Aw + \\Vu\\'^Awt ||(A'u,ut)Aw|| + ||Vm|P ||A?ii||. 


Thus: 


dt 


J-{w, u) 


<C(i?)(l + ||AnJ||). 


Returning to (15.141) . and implementing the above bounds in (15.161) and (15.181) . we have 


{wit),Wt{t)) < C\\{w{to),Wtito))\\Ypi 

ipl 

+ C'(i?,17) [ 

Jto 

< C\\(W{to),Wt{to))\\Ypi 
+ C{R,U) [ 

Jto 


1 + Lf + llwlb + 


k(r)||3dr 


ds 


1 +A:+||w||2 + / \\{w,Wt)\\YpidT 

J s—t* 


ds 


(5.16) 

(5.17) 

(5.18) 


(5.19) 


(5.20) 
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At this point we note by our choice of decomposition—the w problem with null initial data—the initial 
condition = 0 (and hence uJt = 0 as well) on {to — t*,to). Thus: 


[wit),wt{t)) < C\\{w{to),Wt{to))\\Y^i 

ypl 


+ C{R,UX)[^ + K+ sup ||(W,Wt)||y,,] f 

sG [io ^ to 


< c 


[w{to),Wt{to)) 


Ypi 


(5.21) 


(5.22) 


1 _ p-‘^{k,K){t-to) 

+ C{R,U,t*)[l + K+ sup ||(w,uJt)||yp,]--' 

sG[to.t] u}{k,K) 


A global-in-time bound on {wt,wtt) = {w,wt) in Ypi follows by (i) taking supremums in t € [to,c»), then (ii) 
(possibly) up-scaling min[A:, A'] sufficiently (and thus scaling uj{k,K)) so as to absorb the supremum on the 
RHS. This completes the proof of Lemma [5731 


5.5 Decomposition of Flow Dynamics 

Consider data yo €Yp. At a time to sufficiently large, we employ the reduction result Theorem 12.51 and utilize 
initial data xo = {u{to), Ut{to)', u*°) coming from the full flow plate dynamics, and denote the resulting semiflow 
from (|2.10D by Tt{xo)- 

We consider the decomposition of the dynamics Tt{xo) = {u{t),ut{t);u*) into z and w components as in 
(EH) and (15.21) . We then consider the potential flow equation, with w given from (15.2p and initial time to 
sufficiently large. 

r {dt + = Acj)'^ in 

< = -[wt + Uwx]ext in (5.23) 

[((.-(to)=0; (t>nto) = 0. 

In this case, via the smoothing in LemmaEH the Neumann data [wt + Uwxlext G H^{^1) is sufficiently regular 
for Theorem l2.1l to apply. This yields a resulting potential flow {(j)'^,4>f) € C{0, oo; 7L^(]Ri)_) x L 2 (R+))- Owing 
to the calculations in [14] and (9] p. 241] (due to the zero flow initial condition), we have that: 

rntr[[dt + = -{dt + Udx)w - q{w*). (5-24) 

With this identity, we can consider the non-autonomous semiflow S'^{-) on Yp, in reference to the flow-plate 
problem: 

'wtt + kwt + A^w = Po + {dt + Udx)(j)'^\^ - R{w,u) -\- f3z in O 
w = d^w = 0 on r 

^w{to) = t), Wt{to)=0 

' {dt + Udxfcj)'^ = A4>'^ inR3x(0,r), ’ 

r{to) = 0\ </<r(to) = o, 

d^fT = -[{dt + Udx)w\^^^ 

Note that T{w,u) is linear in w, and we have the estimate 

||J-(u;,u)|| < C||Vu||i=(n)||u ;||2 < C\\u\\l\\w\\ 2 . (5.26) 

Owing to the well-posedness of the full flow-plate problem EH [44] and the estimate in (15.261) . one can assert 
well-posedness in Y of the {(j>'^,w) —linear—dynamics by appealing to the results for the linear flow in [HI 
Theorem 2.1, pp. 3128-3128] (see [IS] as well). We now consider the restriction to the flow (j) from the full 
St{yo)- By considering the difference (f> — (j)'^, we can construct the function <f)^, satisfying the corresponding 
potential flow equation with Neumann data 


di/Cj) — {zt UZx\e,yit 
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on K.^. Combining this potential flow and the delayed z plate dynamics, as well as considering the difference 
between the flow-plate systems (ED and (I5.25p . we have constructed a non-autonomous flow /S'f(-) on Yp 
corresponding to the dynamics: 


'ztt + kzt = {dt Udx)4>''\^ - J^iz,u) - j3z in fl 

z = dvZ = 0 on r 

^ z(to) = M(io), ztito) = Ut{to) 

' (dt + Udxfr = in X (0,r), ^ ^ 

(/)^(to) = Hto)-, (io) = 

= - [{dt + Udx)z]^^^ (0’^)- 

Note: by considering k and f3 sufficiently large (as in Lemma [5?^ . we have that (z, z^; z*) decays exponentially 
in Y (uniformity with respect to R). 

Recall that Af corresponds to the points (u, <(>) which satisfy (11.221) . Consider the stationary problem for 
some given function u € HQ(ft): 


{ A^z -h /3z -h R(z,u) = UrQtr[dx4>"'] 
z = d^z = 0 

- U'^dlcj)^ = 0 
dz'f’' = UdxZe^t 

Let 

AT^ = {(z, S HI{Q) X VL^(Ri^) : satisfy 

and let Af^'* C Af^ be the subset of points of the form {0; 0^}. 

5.6 (z, (f>^) Dynamics 

In this section we will show that the desired convergence to equilibria result, as in Theorem 1 1.81 holds for the 
{z,(j)^) dynamics. This will be accomplished by critically using the uniform exponential decay of the (z,Zt) 
plate dynamics, and running an approximation argument with smooth data as in [32]. The main result of this 
section is: 

Theorem 5.6. Assume k > k^ and {3 > fie, with yo G Yp^; then for any ^ -|-oo there is a subsequence 
and a point y^ = (z, 0;^^,0) with (z, 0^) € (i.e., 'z = 0) so that or any p> 0. 

.lim - P'Wy, = 0, 

Proof of Theorem \5. 61 The proof follows through several lemmas. Note that the linear part of (15.271) is equiv¬ 
alent to the system (11.131) (since both have static damping, /3 > 0). The decomposed nonlinear term iF{z,u) 
enjoys the estimate. 


X G n 

xG r 

X G K+ 

X G 9R+ 


(5.28) 


(I5.28P variationally}, 


mz,u)\\ = <C{R,b)\\z\U. (5.29) 

Lemma 5.7. [Uniform-in-Time Hadamard Continuity for S[{-)] Given initial datum = (rt™, uj"; (()™, </))") 
and yo = {uQ,ui;(f>o,(j)i), with corresponding decompositions as in (15.271) and (15.251) . we have that 

sr^p\\Sf{y[[') - S[{yo)\\Y, ^0 as VoUnY 

t>o 

Proof of Lemma [377] With a minor modification to the proof of Theorem [472] in [32] {f{z) replaced by P{z, u)), 
we can obtain global-in-time Hadamard continuity. Indeed, consider the differences z = z"* — z and <i>^ = 
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(j)^ — (j)^ which satisfy the system; 


'ztt + kzt + A^z = {dt + - [T{z"^,u'^) - J^{z,u)] - /3z in 

z = djjZ = 0 on r 

z(to) = ^""(to) - M(io), zt(^o) = u^{to) - Ut{to) 

{dt + in Rl x (0, T), 

$^(to) = (jj’^ito) - (/)(to); ^Uto) = i'Tih) - 
d^(t)^ = - [(ft + t/5x)z]^^^ on X (0,r). 


(5.30) 


We then utilize the velocity multipliers Zt and <i>i (respectively). The linear portion of the dynamics yields the 
standard energy identity for the differences: 

E4t)+ [ k\\ztfdT = E40) + 2um0),z40)) - 2Umt),z4t)) 

Jo 

+ f {E{z^,u"^)-E{z,u),zt)^dT, (5.31) 

Jo 


where we have utilized the notation, 


E.= -{\\Azr 


|ztf+/ 3 ||zf + ||V$f - C/ 2 ||$,f + ||$,f) . 


(5.32) 


By the bounds on Ei^t as in Lemma 11.21 as well as the global-in-time energy bounds iLemma 11.51) . one 
obtains 

Ez{t)+ f fc||zt|pdT 

Jo 

<C{U)EM+C\K{t)f + C [ \{E{z"',un-E{z,u),zt)^\dT. 

Jo 

Via compactness ||zxP < + C'(e)||zp, and hence 

Ez{t)-\- [ k\\zt\\^dr 

Jo 

<C{U)EMEC\\z{t)\\^ + C [ \{J^{z^,u^)-E{z,u),zt)^\dT 

Jo 

Observing that 

< / ||zt||^dr-h ||z(0)|p < / Ez{T)dT + Ez{0), 

Jo Jo 


(5.33) 


ft ^ ft 

Zt(r)dr -I- z( 0 ) 


(5.34) 


we then have: 


E,{t)+ [ k\\ztfdT<C{U)E,{0) + CiR) [ E,{T)dT + C [ \(T{z"^,u"^)-E{z,u),zt}^\dT (5.35) 

Jo Jo Jo 

(where C{R) denotes the dependence of the constant on the size of the ball containing the initial datum ?/° 
and y^). To note finite-in-time Hadamard continuity, let T > 0 be fixed. From (15.311) . we need only consider 
the nonlinear term: ^ 

[ \{E{z'^,u'^) - E{z,u),Zt)\dT. 

Jo 

Note the identity: 

(5- ||Vu'"||2)Az™ - {b-\\Vu\\^)Az = 

6A(z™ -z)- (||Vu'"|p - \\yuf)Az^ - ||Vu|pA(z™ - z) 


31 


(5.36) 









And thus: 


(5.37) 


(5.38) 

(5.39) 


||J-(z™,u'") - T{z,u)\\ < C(6,i?)[||u™ - u \\2 + IN™ - 2 II 2 ]. 

Implementing this bound, we see that 

f {T{z^,u^) - T[z,u),’Zt)dT <C{b,R)f [||u ||2 + ||z|| 2 ]||zt|MT 
^0 Jo 

<Ce{b,R)[[ E^{T)dT+ [ \\u"^-u\\ldT]. 

Jo Jo 

We can now invoke the Hadamard continuity of St{-) for the n-dynamics on any finite [0,T] (Theorem 14.11) . 
as well as the initial conditions in (15.301) . and thus 

\\u^{r)-u{rm<CiT)EM- 

Gronwall’s inequality applied in (I5.35|) with (15.3911 then yields that 

E,{t) < C{b, U, R, T)E,{0). (5.40) 

Hence, for any fixed T, we obtain a Hadamard continuity of the semiflow. 

We now address this continuity on the infinite-time horizon. We return to (15.311) : 

E^{t) + J k\\zt\\^dT < ci^Ez{0)+ \\z{t)\\l+ J - E{z,u),zt)dT 

< c[eM + + \WW + \\E{z^, u^) - T{z, n)||o|| Wllodr} (5.41) 


Again using (I5.34|) . we obtain: 

Ez{t)+J^ k\\ztfdT<c{EM+J^ \\ztfdT + J^ ||J-(z™,M™)-J-(z,n)||o||w||odT} (5.42) 

We invoke the bound (15.3711 and note the global-in-time bounds on the (u™,u("), {u,Ut) and (z’”,z)”), and 
{z, Zt) dynamics in Ypf. 

\\E{z”^,u”^)-E{z,u)\\ < C{b,R)[\\u^ - u\\2 + IN™ - 2II2] < C{R), 


and thus 

Ezit) + J^ k\\ztfdT<c{EM + J^ (INtf+ C'(i?)||zt||)dT} 

Scaling up /3 and k provides exponential decay of zt and z™ to 0, and thus Zt —>■ 0 exponentially—uniform 
with respect to 2/™,2/o G By{R) (see above Lemma [O])- Via this uniform exponential decay, for every e > 0 
there exist a time T* > such tha10 

/ OO 

{W^tf+ C{K)\\zt\\)dT <e. 

Thus, for any t > T* we may write: 

E.it) <c{eM+£ (||zt|N + C(i?)||zt||)dr+^“ {\\ztr + CiR)\\zt\\)dT] (5.43) 

Utilizing (I5.40p . for any e > 0 we have: 

E,it) < C{R,T*){E,{0) + El^^iO)) + e/2 (5.44) 

Taking sufficiently close to will yield that C{R,T*)(^Ez{0) + eI^^{0)) < e/2. This concludes the proof 
of Lemma 15 .71 □ 

^T* depends on the underlying parameters of the problem po, H, U, as well as the time of absorption for and the size of the 
support of the initial flow data p. 
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With uniform-in-time Hadamard continuity available, to continue with an approximation by smooth data 
approach for the dynamics, we must show that the desired convergence to equilibria holds for smooth 

data in the equation. 

Lemma 5.8. Assume that 


xo = (u(to),Mt(to);w*“) e {H^nH^){n) X X L 2 {-t*,0- {H^ n H^){^)) . 

Then the dynamics T^(xo) (such that T“(0) + T^{xo) = Tt{xo), and u = w + z) are uniformly bounded on 
[H'^ n Hq)[Vl) X Hq{Vl) X L 2 (—t*, 0; n Hq){Vl)) . This is to say that: 

Z )\\(^H‘^r\Hg){Q)xHg{Q)xL2{-f‘,0-,{H^nH§){Q)) V t > to- 

Proof of Lemma YSfR We do not repeat the argument in full, as it is directly follows the proof of Theorem 14.71 
mutatis mutandis. Indeed, we consider smooth initial (as above), differentiate the equation (15.11) . and perform 
the Lyapunov analysis in the proof of Theorem (ITTI) . The key fact (which greatly simplifies the proof in this 
case) is below: 


dt 


T{z,u) 


—2(Vu, Vut)Az -I- [6 — llVulpJAzt = 2(Au, ut)Az -I- [6 — ||VM|p]Azt 


(5.45) 


< C(6,i?)[l + ||zi||2]. (5.46) 

0 

Remark 5.2. Note this type of estimate fails for von Karman dynamics because of the precise structure of the 
nonlocality; this is precisely the stage at which our decomposition approach will not obtain for the analogous 
von Karman problem. 

At this point, we repeat the arguments in the proof of the propagation of regularity result Theorem 14.71 
Specifically, the analysis of the term B(t) (involving the decomposition involving terms Pi, Qi) is not needed 
here. In particular, we take Pi and Qi to be zero here. Indeed, utilizing (j5.46|) . we have the estimate: 




—F{z,u),ztt 


^ + C'e||’2tt|P + Ce{R). 


(5.47) 


At the junctures of (I4.17L (I4.28p . and in the proof of Lemma 14.81 we implement (I5.47|) and up-scale k 
appropriately. □ 

Remark 5.3. The subtlety needed in showing the propagation result in Theorem 14.71 linvolving the decomposi¬ 
tion into Pi and Qi) is needed to demonstrate that the propagation of regularity for the original full flow-plate 
dynamics in (u, (f) can be obtained for ANY damping coefficient k > 0. In proving Lemma 15.81 above, we 
already require k large to yield other supporting results, and thus scaling k up to simplify the argument for 
the reduced z dynamics (15.11) analogously is acceptable. 

As a corollary, the global-in-time boundedness in Lemma 15.81 yield the desired convergence to equilibria 
result for smooth initial data to the full 2 problem (15.271) . 

Corollary 5.9. Consider y(f £ ^(T) DY. Thus, for any sequence tn —>• 00 , there is a subsequence depending 
on M (also labeled tn) and a point (z,(j)^) £ (depending on the subsequence, and thus depending on AI) 
such that 

lim I )-i?\\Y„=0, p>0. 

n—foo ^ 

Proof of Corollarv \5.9l Due to the identical structure of the {z, dynamics in ()5.ip to those of (11.11) (with 
f{u) replaced by F(z,u) and removing po): the proof of Theorem 14.31 can be repeated identically for the 
z dynamics dUD, bearing in mind the bound in Lemma 15.81 (Recall: Theorem 14.31 —which translates a 
global-in-time bound on plate dynamics into the desired strong convergence to equilibria for the full flow-plate 
dynamics—does not “care” about the structure of the specific plate equation, only the nature of the coupling. 
See RemarkHH) □ 
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We are now in a position to finish the proof of Theorem 15.61 Now, let e > 0 be given. We consider 
a sequence of smooth initial data G ^(T) with the property that y™ —>■ yo in Y. We may choose M 
sufficiently large so that 

sn^\\S^{y,)-St{y^)\\Y^<e/2. 

t>0 

Hence, there exists a T(e, M) so that for all tn > T{e^ M) 

dy^Sliy^lr) < 

Then 


dy, {SI (yo), r) < dy^ {S^ (yo), (y^)) + dy^ {S^ (y^), r) < 

Owing to the decay of {u,ut) in Lemma lOl we note that, in fact, 
convergence result, and concludes the proof of Theorem 15.61 


e. (5.48) 

This implies the desired 
□ 


5.7 {w,(f)^) Dynamics 

Now we consider the dynamics given by S^{-) on Yp and note the regularity of the dynamics (see Theorem 

lOl) : 

Theorem 5.10. Assume yo G Ypg. For k > kq, the dynamics S'^(-) have the property that for any sequence 
—>■ oo there is a subsequence tn^ and a point y’" = {w, 0; 0) with {w, cjY’) £ Af so that 

.lim \\SZ.{yo) -y^\\y, =0 

for any V p > 0. This implies the desired convergence result for SZ{-) to the set Af 
Proof of Theorem \5.1(A Consider the plate dynamics Tfi-) corresponding to: 

( wtt + (fc + l)wt + A'^w = po + q^{t) - Uwx - \b - ||Vu|p]Aui + fdz in 12 
< u> = dyW = 0 on r (5.49) 

[w(to) = 0, Wt{to)=Q, 10*0=0. 

Lemma [5.3l efuarantees a bound in higher norms H'' for the (lo, Wt', w*) dynamics, satisfying (14.11) in Theorem 
14.31 The analysis of Theorem 14.31 obtains the equivalent of Theorem 14.91 for the w dynamics in (I5.25|) . Here, 
we critically use the fact that the z terms above decay along any sequence of times —>■ oo, i.e., along any 
sequence of times oo we have that z(t„) —>■ 0 in iL^(f2) and Zt{tn) —0 in L 2 {fl). Considering the 

cstirricitjcs 

\\l3z\\o< C||z|| 2 ; [6 - IIVmIPJAw = [6 - ||Vio|p]Ai(; - HVzlpAw, 

and the aforementioned exponential decay, we see that S'^(yo) —t (w, 0; (/>’", 0) (with (io, (/>’") G Af) in Yp for 
any p > 0. This concludes the proof of Theorem 15.101 □ 

5.8 Final Step in the Proof of Theorem 11.81 

Note the fact that for {z, cp^) G Af^’* and {w, 0“) G Af, we have that 

{z,(p^) + {w,(p'^) £ Af (5.50) 

(due the linear structure of the stationary flow equation (11.221) '). By Theorem 15.61 and Theorem 15.101 for any 
sequence —>■ oo, there is a subsequence (also labeled) and two points y^ = (z, 0; (p^,0) (with (S', (p^) £ Af^’*) 
and y’" = {w, 0; 0^", 0) (with {w, ^™) G Af) so that ivo) —t (S', 0; (p^, 0) and (yo) —t {w, 0; 0^", 0) in Yp for 
any p > 0. Thus 

\\Str,{yo) - {P'+ y^)\\y, = \\{SZ{yo) - 1?) + {SZivo) - ir)\\y, (5.51) 

< \\SZiyo) - yZW, + \\SZ,iyo) - y'")\\yp ^ 0, as n-),oo. (5.52) 

Since {z;(p^) + {w](p'^) £ Af, this implies the desired strong convergence result and completes the proof of 
Theorem 11.81 
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6 Open Problem: Eliminating Damping Size Requirements 

In experimental and numerical simulations with panels that have no imposed damping, flutter does not occur 
for subsonic flows (see Section fTTTl) . The techniques herein (and in [T2J [311 [Ml HO]) certainly require some 
damping. However, for finite energy solutions we require a damping coefficient whose size depends on the 
inherent parameters of the problem in order to show convergence of full trajectories to the set of stationary 
states. We note that there is viscous damping in every physical plate system, but we do question whether the 
dependence of the results above on a minimal damping coefficient kmin is necessary. (Especially since there 
is no size restriction for “smooth” initial data—Theorem 14.91 — even for von Karman dynamics.) It seems that 
in a correct model the requirement could be removed in order to yield agreement with experiment, OR that 
there is some physical situation where the structural damping is sufficiently small that the subsonic panel 
may exhibit non-stationary end behavior. The authors could find no discussion of the latter situation in the 
engineering literature. 

In our work above, the key issue necessitating large damping is the interaction of the terms Lu and q{vf) 
with the multiplier analysis in Section 15.31 and Section 15.41 The following intermediate assertion is presently 
clear: Taking b = 0^ po = 0, U << 1, will not eliminate the need for a damping restriction to obtain strong 
stabilization of the full flow-plate dynamics due to the non-conservative, non-dissipative contributions from 
q. An interesting question is: under what restrictions (of the form 6 = 0 and po = Oj or at least small in 
some sense) can we expect stability of the plate dynamics—no non-trivial stationary states? From what is 
known physically, large damping, diminished loading, or small unperturbed flow velocities should lead to a 
simplification of the stationary set Af. 
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